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ABSTRACT 



We report on the 7-ray activity of the blazar Mrk 501 during the first 480 days 
of Fermi operation. We find that the average LAT 7-ray spectrum of Mrk 501 
can be well described by a single power-law function with a photon index of 
1.78 ±0.03. While we observe relatively mild flux variations with the Fermi-LAT 
(within less than a factor of 2), we detect remarkable spectral variability where 
the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, 
and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not 
correlate with the measured flux variations above 0.3 GeV. In this paper, we also 
present the first results from the 4.5-month-long multifrequency campaign (2009 
March 15 - August 1) on Mrk 501, which included the VLBA, Swift, RXTE, 
MAGIC and VERITAS, the F-GAMMA, GASP-WEBT, and other collabora- 
tions and instruments which provided excellent time and energy coverage of the 
source throughout the entire campaign. The extensive radio to TeV data set from 
this campaign provides us with the most detailed spectral energy distribution yet 
collected for this source during its relatively low activity. The average spectral 
energy distribution of Mrk 501 is well described by the standard one- zone syn- 
chrotron self-Compton model. In the framework of this model, we find that the 
dominant emission region is characterized by a size < 0.1 pc (comparable within 
a factor of few to the size of the partially- resolved VLBA core at 15-43 GHz), and 
that the total jet power (~ 10 44 ergs _1 ) constitutes only a small fraction (~ 10~ 3 ) 
of the Eddington luminosity. The energy distribution of the freshly-accelerated 
radiating electrons required to fit the time-averaged data has a broken power-law 
form in the energy range 0.3 GeV— 10 TeV, with spectral indices 2.2 and 2.7 below 
and above the break energy of 20 GeV. We argue that such a form is consistent 
with a scenario in which the bulk of the energy dissipation within the domi- 
nant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We 
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find that the ultrarelativistic electrons and mildly relativistic protons within the 
blazar zone, if comparable in number, are in approximate energy equipartition, 
with their energy dominating the jet magnetic field energy by about two orders 
of magnitude. 

Subject headings: acceleration of particles — radiation mechanisms: non-thermal 
- galaxies: active — BL Lacertae objects: general — BL Lacertae objects: 
individual (Mrk501) — gamma rays: observations — radio continuum: galaxies 
— ultraviolet: galaxies — X-rays: galaxies 



1. Introduction 



Blazars constitute a subclass of radio- loud active galactic nuclei (AGN), in which a 
jet of magnetized plasma assumed to emanate with relativistic bulk velocity from close to 
a central supermassive black hole points almost along the line of sight. The broadband 
emission spectra of these objects are dominated by non-thermal, strongly Doppler-boosted 
and variable radiation produced in the innermost part of the jet. Most of the identified 
extragalactic 7-ray sources detected with the EGRET instrument (Hartman et al. 1999) on 
board the Compton Gamma Ray Observatory belong to this category. Blazars include flat- 
spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs). Even though blazars 
have been observed for several decades at different frequencies, the existing experimental 
data did not permit unambiguous identification of the physical mechanisms responsible for 
the production of their high-energy (7-ray) emission. Given the existing high- sensitivity 
detectors which allow detailed study of the low-energy (synchrotron) component of blazar 
sources (extending from radio up to hard X-rays), one of the reasons for the incomplete 
understanding of those objects was only moderate sensitivity of previous 7-ray instruments. 
This often precluded detailed cross-correlation studies between the low- and high-energy 
emission and did not provide enough constraints on the parameters of the theoretical models. 
Some of the open and fundamental questions regarding blazar sources are (i) the content of 
their jets, (ii) the location and structure of their dominant emission zones, (iii) the origin of 
their variability, observed on timescales from minutes to tens of years, (iv) the role of external 
photon fields (including the extragalactic background light, EBL) in shaping their observed 
7-ray spectra, and (v) the energy distribution and the dominant acceleration mechanism for 
the underlying radiating particles. 



The Large Area Telescope (LAT) instrument (lAtwood et al.l 120091 ) on board the Fermi 
Gamma-ray Space Telescope satellite provides a large improvement in the experimental ca- 
pability for performing 7-ray astronomy, and hence it is shedding new light on the blazar 
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phenomenon. In this paper, we report on the Fermi observations of the TeV-emittin g high- 
frequency-peaked — or, according to a more recent classification (lAbdo et al.ll2010d j. high- 
synchrotron-peaked (HSP) — BL Lac object Markarian501 (Mrk501; RA=16 /l 45 m 52.22 s , 
Dec= 39° 45' 36.6" , J2000, redshift z = 0.034), which is one of the brightest extragalactic 
sources in the X-ray/TeV sky. Mrk 501 was the second extragalactic ob ject (after Markar 



ian421) identified as a very high energy (thereafter VHE) 7-ray emitter (IQuinn et al.l 11996 



Bradbury et al.lll997l ). After a phase of moderate emission lasting for about a year following 



its discovery (1996), Mrk 501 went into a state of surprisingly high activity and strong vari 
ability, becoming > 10 times br ighter than the Crab Nebula at energies > 1 TeV, as reported 
by var ious instruments /groups ( Catanese et a0l997 ; Samuelson et al. 1998 ; Aharonian et al. 



1999aH 3: iDjannati-Atai et al 



199% . In 1998-1999, the mean VHE 7-ray flux dropped by 



an order of magnitud e, and the overall VHE spectrum softened significantly (jPironl 12000 



Aharonian et al. 



20011 ). In 2005, 7-ray flux variability on minute timescales was observed in 



the VHE band, thus establis hing Mrk 501 as one of the sources with the fastest 7-ray flux 
changes (lAlbert et al.ll2007al ). During the 2005 VHE flux variations (when Mrk 501 was 3-4 
times dimmer than it was in 1997), significant spectral variability was detected as well, with a 
clear "harder when brighter" behavior. Those spectral variations are even more pronounced 
w hen compared wi t h the spectrum measured during the low activity level recently reported 
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, but the nature of this object 



is still far from being understood. The main reasons for this lack of knowledge are the sparse 
multifrequency data during long periods of time, and the moderate sensitivity available in 
the past to study the 7-ray emission of this source. Besides, most of the previous multifre- 
quency campaigns were triggered by an enhanced flux level in some energy band, and hence 
much of our information about the source is biased towards "high-activity" states, where 
perhaps distinct physical processes play a dominant role. In addition, until now we knew 
very little about the GeV emission of Mrk 501. The only detection reported at GeV energies 
before Fermi was in lKataoka et al.l (Il999l ). but t he significance of this detection was too low 
to include Mrk 501 in the 3rd EGRET catalog jHartman et al.l Il999h . Moreover, Mrk 501 
was not detected by E GRET during the large X-ray and VHE 7-ray flare which lasted for 
several months in 1997 (IPian et al.lll998l ). 

The large improvement in the performance provided by the Fermi-LAT compared with 
its predecessor, EGRET, provides us with a new perspective for the study of blazars like 
Mrk 501. However, it is important to emphasize that blazars can vary their emitted power 
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by one or two orders of magnitude on short timescales, and that they emit radiation over 
the entire observable electromagnetic spectrum (from ~ 10 _6 eV up to ~ 10 13 eV). For this 
reason, the information from Fermi-LAT alone is not enough to understand the broadband 
emission of Mrk501, and hence simultaneous data in other frequency ranges are required. In 
particular, the frequency ranges where the low- and high-energy spectral components peak 
in the SED representation are of major importance. In the case of Mrk501, those peaks 
are typically located around IkeV (low-energy bump) and 100 GeV (high-energy bump), 
and hence simultaneous UV/X-ray and GeV/TeV observations are essential for the proper 
reconstruction of the overall SED of Mrk501. At TeV energies there has been a substantial 
improvement in the instrumental capability as a result of the deployment of a new gener- 
ation of imaging atmospheric Cherenkov telescopes (IACTs). In particular, for the study 
of Mrk501, the new telescope systems MAGIC and VERITAS provide greater sensitivity, 
wider energy range and improved energy resolution compared with the previous generation of 
instruments. Simultaneous observations with Fermi-LAT and IACTs like MAGIC or VER- 
ITAS (potentially covering six decades in energy, from 20 MeV to 20 TeV) can, for the first 
time, significantly resolve both the rising and the falling segments of the high-energy emis- 
sion component of Mrk501, with the expected location of the SED peak in the overlapping 
energy range between those instruments. Because of the smaller collection area, and the 
self-veto proble the sensitivity of EGRET to detect 7-rays with energies larger than 10 
GeV was about two orders of magnitude lower than that of Ferrm-LAT^. Besides, during the 
period of operation of EGRET, the sensitivity of the previous generation of IACTs was only 
moderate, with relatively low sensitivity below 0.5 TeV. Therefore, the higher sensitivity and 
larger energy range of the newer 7-ray instruments has become a crucial tool for studying 
Mrk501, and the blazar phenomenon in general. 

In order to exploit the performance of the Fermi-LAT and the new IACTs, as well 
as the capabilities of several existing instruments observing at radio-to-X-ray frequencies, 
a multifrequency (from radio to TeV photon energies) campaign was organized to monitor 
Mrk501 during a period of 4.5 months, from mid-March to August 2009. The scientific 
goal was to collect a very complete, simultaneous, multifrequency data set that would allow 
current theoretical models of broadband blazar emission to be tested. This, in turn, should 
help us to understand the origin of high-energy emission of blazar sources and the physi- 
cal mechanisms responsible for the acceleration of radiating particles in relativistic jets in 



1 The self- veto problem in EGRET is the degradation of the effective area at high energies (>5 GeV) due 
to backsplash of secondary particles from the calorimeter causing the anticoincidence system to veto the 
event. This problem is substantially reduced in LAT by using a segmented anticoincidence detector. 

2 This estimate includes the larger exposure from Fermi-LAT due to the 4 times larger field of view. 
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general. In this paper, the only reported result from the multifrequency observations is the 
overall SED averaged over the duration of the observing campaign. A more in-depth analysis 
of the multifrequency data set will be given in a forthcoming paper. The scientific results 
from the data collected during the two-day time interval 2009 March 23-25 (which includes 
extensive observation s with the Suzaku X-ray satellite) will be reported in a separate paper 
( lAcciari et al.ll2010bl ). The paper is organized as follows. In §2 we introduce the LAT instru- 
ment and describe the LAT data analysis. In §3 we report on the flux/spectral variability of 
Mrk501 observed during the first 16 months of Fermi-LAT operation, a nd compare it with 
the flux variability observed in X-rays by the all-sky instru ments RXTE (IBradt et al.lll993f ) 
All Sky Monitor (ASM) and the Swift f lGehrels et al1l2004j ) Burst Alert Telescope (BAT). In 
§4 we analyze the 7-ray spectrum of Mrk501 measured by Fermi-LAT in the energy range 
0.1 — 400 GeV. §5 reports on the overall SED obtained during the 4.5-month-long multifre- 
quency campaign organized in 2009. §6 is devoted to SED modeling, the results of which 
are further discussed in §7. Conclusions are presented in §8. 



2. Fermi-LAT Data Selection and Analysis 

The Fermi-LAT is an instrument to perform 7-ray astronomy above 20 MeV. The in- 
strument is an array of 4 x 4 identical towers, each one consisting of a tracker (where the 
photons are pair-converted) and a calorimeter (where the energies of the pair-converted pho- 
tons are measured). The entire instrument is covered with an anticoincidence detector to 
reject charged-particle background. The LAT has a peak effective area of 0.8 m 2 for 1 GeV 
photons, an energy resolution typically better than 10% and a field of view (FoV) of about 
2.4 sr, with an angular resolution (68% containment angle) better t han 1 ° for energies above 



1 GeV. Further details on the LAT can be found in lAtwood et al.l (120091 ) . 



The LAT data reported in this paper were collected from 2008 August 5 (MJD 54683) 
to 2009 November 27 (MJD 55162). During this time, the Fermi-LAT instrument operated 
mostly in survey mode. The analysis was performed with the Fermi Science Tools software 
package version v9rl5p6. Only events with the highest probability of being photons - 
those in the "diffuse" class — were used. The LAT data were extracted from a circular 
region of 10° radius centered at the location of Mrk501. The spectral fits were performed 
using photon energies in the energy range 0.3 — 400 GeV. At photon energies above 0.3 GeV 
the effective area of the instrument is relatively large (> 0.5 m 2 ) and the angular resolution 
relatively good (68% containment angle smaller than 2°). In particular, because of the better 
angular resolution, the spectral fits using energies above 0.3 GeV (instead of 0.1 GeV) are 
less sensitive to possible contamination from unaccounted (perhaps transient), neighboring 
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7-ray sources and hence have smaller systematic errors, at the expense of reducing somewhat 
the number of photons from the source. In addition, a cut on zenith angle (> 105°) was 
applied to reduce contamination from Earth-albedo 7-rays, which are produced by cosmic 
rays interacting with the upper atmosphere. 

The background model used to extract the 7-ray signal includes a Galactic diffuse emis- 
sion component and an isotropic component. The model that we adopted for the Galactic 
component is gll_iem_ v02.fit§ The isotropic component, which is the sum of the extra- 
galactic diffuse emission and the residual charged-particle background, is parametrized here 
with a single power-law function. To reduce systematic uncertainties in the analysis, the 
photon index of the isotropic component and the normalization of both components in the 
background model were allowed to vary freely during the spectral point fitting. Owing to 
the relatively small size of the region analyzed (radius 10°) and the hardness of the spec- 
trum of Mrk501, the high-energy structure in the standard tabulated isotropic background 
spectrum isotropic_iem_v02.txt does not dominate the total counts at high energies. In ad- 
dition we find that for this region a power-law approximation to the isotropic background 
results in somewhat smaller residuals for the overall model, possibly because the isotropic 
term, with a free spectral index, compensates for an inaccuracy in the model for the Galac- 
tic diffuse emission, which is also approximately isotropic at the high Galactic latitude of 
Mrk501 (b ~ 39°). In any case, the resulting spectral fits for Mrk501 are not significantly 
different if isotropic_iem_v02.txt is used for the analysis. In additi on, the model also in- 



cludes five nearby sources from the 1FGL catalog flAbdo et alJl2010b ): 1FGL J1724.0+4002, 
1FGL J1642. 5+3947, 1FGL J1635. 0+3808, 1FGL J1734. 4+3859, and 1FGL J1709. 6+4320. 
The spectra of those sources were also parameterized by a power-law functions, whose pho- 
ton index values were fixed to the values from the 1FGL catalog, and only the normaliza- 
tion factors for the single sources were left as free parameters. The spectral analysis was 
performed with the post-launch instr ument-response func tions P6_V3_DIFFUSE using an un- 



binned maximum-likelihood method (jMattox et al.lll996l ). The systematic uncertainties on 
the flux were estimated as 10% at 0.1 GeV, 5% at 560 MeV and 20% at lOGeV and abovr^ 



3. Flux and Spectral Variability 

The high sensitivity and survey-mode operation of Fermi-LAT permit systematic, un- 
interrupted monitoring of Mrk501 in 7-rays, regardless of the activity level of the source. 



• fottp : //f ermi . gsf c .n asa. gov/ssc/data/access/lat/BackgroundModels .html 



http : //f ermi . gsf c .nasa. gov/ssc/data/analysis/LAT_caveats .html 
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The measured 7-ray flux above 0.3 GeV and the photon index from a power-law fit are 
shown in the left panel of Figured] The data spans the time from 2008 August 5 (MJD 
54683) to 2009 November 27 (MJD 55162), binned in time intervals of 30 days. The Test 
Statistic (TS) valuel for the 16 time intervals are all in excess of 50 (i.e., ~ 7 standard 
deviations, hereafter cr), with three-quarters of them greater than 100 (i.e., ~ 10 cr). During 
this 480-day period, Mrk501 did not show any outstanding flaring activity in the Fermi-LAT 
energy range, but there appear to be flux and spectral variations on timescales of the order 
of 30 days. During the 120-day period MJD 54862-54982, the photon flux above 0.3 GeV 
was (3.41 ± 0.28) x 10 -8 phcm~ 2 s , which is about twice as large as the averaged flux 
values before and after that time period, which are (1.65 ± 0.16) x 10~ 8 phcm -2 s _1 and 
(1.84 ± 0.17) x 10~ 8 phcm -2 s -1 , respectively. Remarkably, the photon index changed from 
2.51 ± 0.20 for the first 30-day interval of this "enhanced-flux period" to 1.63 ± 0.09 for 
the last 30-day interval. As shown in the red legend of the bottom plot in the left panel 
of Figure[[J a constant fit to the photon index values of this 120-day period gives a null 
probability of 10 -4 , hence a deviation of 4 a. A constant fit to the entire 16-month period 
gives a null probability of 2.6 x 10 -3 , hence spectral variability is detected for the entire 
data set at the level of 3 cr. It is worth stressing that the spectral variability in the 480-day 
time interval is entirely dominated by the spectral variability occurring during the 120-day 
time interval of MJD 54862-54982, with no significant spectral variability before or after 
this "enhanced-flux period" . The right plot in Figured] does not show any clear correlation 
between the flux and the spectral vari ations. The discrete correlation function computed as 
prescribed in lEdelson &: Kroliki (119881 ) gives DCF = 0.5 ± 0.3 for a time lag of zero. 



Mrk501 is known for showing spectral variability at VHE 7-ray energies. During the 
large X-ray /7-ray flare in 1997, Whi pple and (especially) CAT observations showed ev idence 
of spectral curvature and variability (ISamuelson et al.lll998l ; iDjannati-Ata'i et al.lll999l ). The 
spectral changes are larger when comparing the measure ments from 1997 with the low state s 
from 1998 and 1999, as reported by CAT and HEGRA JPironlbood : lAharonian et al.ll200lh . 
The MAGIC telescope, with lower energy threshold and higher sensitivity than the Whipple, 
HEGRA and CAT telescopes, observed remarkable spectral variability in 2005, when the 7- 



ray ac tivity of Mrk501 was significantly smaller than that observed in 1997 ([Albert et al. 
2007a|). The spectral variability is even larger when comparing the MAGI C measurements 



from 2005 with those from 2006 when the source was in an even lower state (Anderhub et al. 



20091 ). However, despite the measured spectral variability at VHE 7-ray energies, the out- 



5 The Test Statistic value quantifies the probability of having a point 7-ray source at the location specified. 
It is roughly the squar e of the significance value: a TS of 25 would correspond to a signal of approximately 
5 standard deviations (jMattox et al.lll996l ). 
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standing spectral steepening at GeV energies observed during the time interval MJD 54862- 
54892 was not envisioned in any of the previous works in the literature; the modeled spectrum 
of Mrk501 at GeV energies was always assumed to be hard (photon indices ~ 1.5—1.8). This 
observational finding, further discussed in §Hand §0 shows the importance of having a 7-ray 
instrument capable of long-term, uninterrupted, high-sensitivity monitoring of Mrk501 and 
other HSP BL Lacs, and it points to the important role Fermi-LAT will play in improving 
our understanding of the physics behind the blazar phenomenon. 

The Fermi-LAT capability for continuous source monitoring is complemented at X-ray 
frequencies by RXTE-ASM and Swift-BAT , the two all-sky instruments that can probe the 
X-ray activity of Mrk501 on a 30-day timescale. Figure[2] shows the fluxes measured by 
ASM in the energy range 2 — lOkeV, by BAT in the energy range 15 — 50keV, and by 
LAT in two different energy bands: 0.2 — 2 GeV (low-energy band) and > 2 GeV (high- 
energy band)B The data from RXTE-ASM were obtained from the ASM web pag fl The 
data were filtered according to the prescription provided there, and the weighted average 
over all of the dwella^l was determined for the 30-day time intervals defined for the Fermi 
data. The data from Swift-BAT were gathered from the BAT web page^l We retrieved the 
daily averaged BAT values and made the weighted average over all the days from the 30- 
day time intervals defined for the Fermi data. The X-ray intensity from Mrk501, averaged 
over the 16 months, is 0.25 ± 0.01 cts^ 1 per Scanning Shadow Camera (SSC) in ASM, and 
(0.52 ± 0.05) x lO^cts^cm- 2 in BAT (close to the BAT 30-day detection limit). This 
X-ray activity is compatible with that recorded in recent years, but quite different from the 
activity of the source during 1997, when the ASM flux was above lets -1 per SSC during 
most of the year, with a peak well above 2 ct s _1 around June 1997. 

As noted previously (51]), Mrk501 is not in the 3rd EGRET catalog, although there was 
a marginally significa nt EGRET detection during the 7-ray outburst (with no clear X-ray 



counterpart) in 1996 (IKataoka et al.lll999l ). At that time, the source was detected at a level 
of 4.0 a at energies above 0.1 GeV and at 5.2 o above 0.5 GeV. The flux from the EGRET 
1996 flare above 0.5 GeV was (6 ±2) x 10~ 8 phcm -2 s -1 , which is about five times higher than 
the average flux observed by Fermi from 2008 August 5 (MJD 54683) to 2009 November 27 



6 The fluxes depicted in the Fermi-LAT light curves were computed fixing the photon index to 1.78 
(average index during the first 480 days of Fermi operation) and fitting only the normalization factor of the 
power-law function. 

7 http : //xte .mit . edu/ASM_lc . html 

8 A dwell is a scan/rotation of the ASM Scanning Shadow Cameras lasting 90 seconds. 
s http : //swift .gsfc.nasa.gov/docs/swift/results/transients/ 
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(MJD 55162), namely (1.39 ±0.07) x 10~ 8 ph cm" 2 s" 1 (also above photon energy 0.5 GeV). 
The Fermi-LAT flux measured during the 120 days with the "enhanced" 7-ray activity (MJD 
54862-54982) is (2.03 ± 0.18) x lO^phcm- 2 s" 1 (above photon energy 0.5 GeV), about a 
factor of three lower than that detected by EGRET in 1996. 

In spite of the relatively low activity, the ASM and BAT fluxes show some flux variations 
and a positive correlation between the fluxes measured by these two instruments. The 
discrete correlation function for the ASM/BAT data points shown in Figure^] is DCF = 
0.73 ± 0.17 for a time lag of zero. On the other hand, the X-ray ASM/BAT fluxes are 
not significantly correlated with the 7-ray LAT fluxes. We found, for a time lag of zero, 
DCF = 0.32 ±0.22 for the ASM/LAT (< 2 GeV) and DCF = 0.43 ±0.30 for the ASM/LAT 
(> 2 GeV) flux data points shown in Figure[2J It is also interesting to note that the largest 
flux variations occur at the highest Fermi energies (> 2 GeV), where the 7-ray flux increased 
by one order of magnitude during the 120-day interval MJD 54862-54892. This trend is 
consistent with the photon index hardening revealed by the spectral analysis reported above 
(see Figure[T]). 



We followed the description given in lVaughan et all (120031 ) to quantify the flux variability 



by means of the fractional variability parameter, F var , as a function of photon energy. In 
order to account for the individual flux measurement e rrors (a P . r r,i), we used the "excess 



varia nce" as an estimator of the intrinsic source variance ( INandra et ajJ[l997J ; lEdelson et al. 



20021 ). This is the variance after subtracting the expected contribution from the measurement 



errors. For a given energy range, F var is calculated as 



S 2 (cTg rr ) 



2 



where (F 7 ) is the mean photon flux, S is the standard deviation of the N flux points, and 
(<Tg rr ) is the average mean square error, all determined for a given energy bin. 

Figure|3] shows the F var values derived for the four different energy ranges and the time 
window covered by the light curves shown in Figure|2j The source is variable at all energies. 
The uncertainty in the variability quantification for the Swift-BAT energies is large due to 
the fact that Mrk501 is a relatively weak X-ray source, and is therefore difficult to detect 
above 15keV in exposure times as short as 30 days. On the contrary, the variability at the 
RXTE-ASM and, especially, Fermi-LAT energies, is significant (> 3 a level). The amplitude 
variability in the two X-ray bands is compatible within errors, and the same holds for the 
variability in the two 7-ray bands. As shown in Figure[3j for the hypothesis of a constant 
F var over the four energy bands one obtains x 2 = 3.5 for 3 degrees of freedom (probability 
of 0.32), implying that the energy- dependent variability is not statistically significant. It 
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is worth noticing that the limited sensitivity of ASM and (particularly) BAT instruments 
to detect Mrk501 in 30-day time intervals, as well as the relatively stable X-ray emission 
of Mrk501 during the analyzed observations, precludes any detailed X-ray/7-ray variability 
and correlation analysis. 



4. Spectral Analysis up to 400 GeV 

The large effective area of the Fermi-LAT instrument permits photon energy recon- 
struction over many orders of magnitude. As a result, the spectrum of Mrk501 could be 
resolved within the energy range 0.1 — 400 GeV, as shown in FigureHl This is the first 
time the spectrum of Mrk501 has been studied with high accuracy over this energy range. 
The fluxes were computed using the analysis procedures described in £j2j The black line in 
FigureH] is the result of an unbinned likelihood fit with a single power-law function in the 
energy range 0.3 — 400 GeVcj, and the red contour is the 68% uncertainty of the fit. The 
data are consistent with a pure power-law function with a photon index of 1.78 ± 0.03. The 
black data points result from the analysis in differential energy range @ (log = 0.4). 
The points are well within 1 — 2a from the fit to the overall spectrum (black line), which 
confirms that the entire Fermi spectrum is consistent with a pure power-law function. Note, 
however, that, due to the low photon count, the error bars for the highest energy data points 
are rather large. The predicted (by the model for Mrk501) number of photons detected by 
LAT in the energy bins 60 — 160 GeV and 160 — 400 GeV are only 11 and 3, respectively. 
Therefore, even though the signal significance in the highest-energy bins are very high due 
to the very low background (the TS values for the two highest-energy ranges is 162 and 61, 
respectively), the large statistical uncertainties could hide a potential turnover in the spec- 
trum of Mrk501 around 100 GeV photon energies. As we know from past observations, the 



VHE spectrum is significan tly softer than the one observed by Fermi (e.g., lAharonian et al. 



200 ll ; lAnderhub et all 120091 ). and hence the spectrum of Mrk501 must have a break around 



the highest Fermi-LAT energies. 

In §|3]we reported remarkable spectral variability during the 120-day time interval MJD 



10 The unbinned likelihood fit was performed on photon energies above 0.3 GeV in order to reduce system- 
atics. See $2]for further details. 

n Because the analysis was carried out in small energy ranges, it was decided to fix the spectral index 
at 1.78 (the value obtained from fitting the entire energy range) and fit only the normalization factor. We 
repeated the same procedure fixing the photon indices to 1.5 and 2.0 and found no significant change. 
Therefore, the results from the differential energy analysis are not sensitive to the photon index used in the 
analysis. 
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54862-54982, when Mrk501 was characterized by a photon flux (at > 0.3 GeV) twice as 
large as during the rest of the exposure. In order to understand better the behaviour of the 
source during that time, we produced SED plots (analogous to that of FigureS]) for each of 
the 30-day time intervals from the period with the enhanced flux level. These are shown in 
Figure[5j together with the SED plots from the 30-day time intervals before and after this 
120-day epoch, which are representative of the average source behaviour during the other 
360 days. The variability of the SED data points below a few GeV is rather mild (factor of 
two), but above a few GeV the spectra vary substantially (factor of ten). The 7-ray signal 
at the highest energies is suppressed during MJD 54862-54982, while it increases by a large 
factor during MJD 54952-54982, where the analysis model for Mrk501 predicts 2.0 photons 
in the energy range 160 — 400 GeV. It is worth stressing that for the SED from FigurelU 
which corresponds to the total exposure of 480 days, the analysis model for Mrk501 predicts 
only 3.2 photons in the highest energy bin. Hence the time interval MJD 54952-54982 
holds almost all the signal detected by LAT in the energy range 160 — 400 GeV during 16 
months. The situation changes somewhat for the lower energy bin 60 — 160 GeV, for which 
the analysis model for Mrk501 predicts 2.4 photons for the time interval MJD 54952-54982, 
while it does predict 11.3 photons for the entire 16-month time interval. Fortunately, the 
30-day time interval characterized by hard spectrum is covered by the 4. 5- mo nth campaign 
that we organized, and hence simultaneous multifrequency observations (radio to TeV) are 
available for this particular period, as discussed further below. 



5. Broadband Spectral Energy Distribution of Mrk501 

As mentioned in £TIJ we organized a multifrequency campaign (from radio to TeV photon 
energies) to monitor Mrk501 during a time period of 4.5 months. The observing campaign 
started on March 15, 2009 (MJD 54905) and finished on August 01, 2009 (MJD 55044). The 
observing goal for this campaign was to sample the broadband emission of Mrk501 every 
5 days, which was largely accomplished whenever the weather and/or technical limitations 
allowed. The underlying scientific goal has already been outlined in §1. A detailed analysis of 
the multifrequency variability and correlations, as well as the evolution of the overall spectral 
energy distribution with time, will be reported in a forthcoming paper. In this section of the 
manuscript, we describe the source coverage during the campaign and the data analysis for 
several of the participating instruments, and we report on the averaged SED resulting from 
the campaign. The modeling of these data and the physical implications are given in §6] and 
below, respectively. 
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5.1. Details of the Campaign: Participating Instruments and Temporal 

Coverage 

The list of all the instruments that participated in the campaign is given in Tabled] and 
the scheduled observations can be found onlinj^l. In some cases the planned observations 
could not be performed due to bad observing conditions, while in some other occasions the 
observations were performed but the data could not be properly analyzed due to technical 
problems or rapidly changing weather conditions. In order to quantify the actual time and 
energy coverage during the campaign on Mrk501, Figure[6] shows the exposure time as a 
function of the energy range for the instruments/observations used to produce the SED 
shown in Figure[SJ Apart from the unprecedented energy coverage (including, for the first 
time, the GeV energy range from Fermi-LAT) , the source was sampled quite uniformly with 
the various instruments participating in the campaign and, consequently, it is reasonable to 
consider the SED constructed below as the actual average (typical) SED of Mrk501 during 
the time interval covered by this multifrequency campaign. The largest non-uniformity in the 
sampling of the source comes from the Cherenkov Telescopes, which are the instruments most 
sensitive to weather conditions. Moreover, while there are many radio/optical instruments 
spread all over the globe, there are only three Cherenkov Telescope observatories in the 
northern hemisphere we could utilize (MAGIC, VERITAS, Whipple). Hence, the impact of 
observing conditions was more important to the coverage at the VHE 7-ray energies. 

We note that FigureE] shows the MAGIC, VERITAS and Whipple coverage at VHE 
7-ray energies, but only the MAGIC and VERITAS observations were used to produce 
the spectra shown in FigureEJ The more extensive (120 hr), but less sensitive, Whipple 
data (shown as gr ey boxes in Figure^]) were primarily taken to determine the light curve 



( IPichel et al.ll2009l ) and a re-optimization was required to derive the spectrum which will be 



reported elsewhere. 

In the following paragraphs we briefly discuss the procedures used in the data analysis 
of the instruments participating in the campaign. The analysis of the Fermi-LAT data was 
described in §2] and the results obtained will be described in detail in §5.21 



5.1.1. Radio Instruments 



Radio data were taken for this campaign from single-dish telescopes, one mm-interferometer, 
and one Very Long Baseline Interferometry (VLBI) array, at frequencies between 2.6 GHz 



12 tittp s : //confluence . slac . Stanford. edu/display/GLAMCDG/Campaign+on+Mrk501+(March+2009+ to+July+ 2009) 
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and 225 GHz (see Tabled]). The single-dish telescopes were the Effelsberg 100 m radio tele- 
scope, the 32 m Medicina radio telescope, the 14 m Metsahovi radio telescope, the 32 m Noto 
radio telescope, the Owens Valley Radio Observatory (OVRO) 40 m telescope, the 26 m 
University of Michigan Radio Astronomy Observatory (UMRAO) and the 600 meter ring 
radio telescope RATAN-600. The mm-interferometer is the Sub-millimeter Array (SMA). 
The NRAO Very Long Baseline Array (VLBA) was used for the VLBI observations. For 
the single-dish instruments and SMA, Mrk501 is point-like and unresolved at all observing 
frequencies. Consequently, the single-dish measurements denote the total flux density of the 
source integrated over the wh ole source extension. D e tails of the obse r ving strategy and 
data red uction can be found inlFuhrmann et all (120081); lAngelakis et al.l (120081 F- GAMMA 



projec t), iTerasranta et all (Il998l Metsahovi), 
( boOll Medicina and Noto). lKovalev et all (11999 
ration, OVRO. 



Aller et all (119851 . UMRAO), 



RATAN-600) and lRichards et al 



Venturi et al 



m prepa- 



In the case of the VLBA, the data were obtained at various frequencies from 5 GHz 
to 43 GHz through various programs (BP143, BK150 and MOJAVE). The data were re- 
duced following s tandard procedures for data reduction and calibration (see, for example, 
Lister et all 120091 for a description of the MOJAVE program which provided the 15 GHz 
data). Since the VLBA angular resolution is smaller than the radio source extension, mea- 
surements were performed for the most compact core region, as well as for the total radio 
structure at parsec scales. The VLBA core size was determined with two-dimensional cir- 
cular or elliptical Gaussian fits to the measured visibilities. The FWHM size of the core 
was estimeted to be in the range 0.14-0.18 mas at the highest observing frqeuencies, 15-43 
GHz. Both the total and the core radio flux densities from the VLBA data are depicted in 
Figure[HJ 



5.1.2. Optical and Near-IR Instruments 



The coverage at optical frequencies was obtained through various telescopes around 
the globe, and this decreased the sensitivity to weather /technical difficulties and provided 
good overall coverage of the source, as depicted in F igureED Many of th e obse rvations were 
performed within the GASP-WEBT program (e.g., Ivillata et aD I2008I . l2009h : that is the 
case for the data collected by the telescopes at Abastumani, Lulin, Roque de los Muchachos 
(KVA), St. Petersburg, Talmassons, and Valle d'Aosta observatories (R band), and also 
for Campo Imperatore (near- infrared frequencies, JHK bands). In addition, the telescopes 
GRT, ROVOR and MitSume provided data with various optical filters, while OAGH and 
WIRO provided data at near- infrared wavelengths. See Table[T] for further details. 
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All the instruments used the calibration stars reported in IVillata et al.l (ll998[). an d the 
Galactic extinction was corrected with the coefficients given in ISchlegel et al.l (ll998f ). On 
the other hand, the flux from the host gala xy, which in the R band accounts for about two- 
thirds of the overall measured optical flux (INilsson et al.ll2007| ). was not subtracted. As can 
be seen from FigureEJ the host galaxy contribution shows up as an additional (narrow) bump 
in the SED with the peak located at infrared frequencies and the flux decreasing rapidly with 
increasing frequency. At frequencies above 10 15 Hz, the blazar emission again dominates the 
radiative output of Mrk501. 



5.1.3. Swift-UVOT 



The Swift-XJltra,- Violet/Optical Telescope (UVOT; iRoming et al.l 120051 ) data used in 
this analysis include all the observations performed during the time interval MJD 54905 
and 55044, which amounts to 41 single pointing observations that were requested to provide 
UV coverage during the Mrk501 multifrequency campaign. The UVOT telescope cycled 
through each of six optical and ultraviolet passbands (V, B, U, UVW1, UVM2, UVW2). 
Photometry was computed using a 5arcsec source region around Mrk501 usi ng a custom 



UVO T pipeline that obtains similar photometric results to the public pipeline ( IPoole et al. 



20081 ). The custom pipeline also allows for separate, observation-by-observation corrections 
for astrometric mis- alignments ( Acciari et al.ll2010d . in preparation). A visual inspection 



was also performed on each of the observations to ensure proper data quality selection and 
correction. The flu x measurements ob tained have been corre cted for Galactic extinction 
E B _ V = 0.019 mag (ISchlegel et al.lll998l ) in each spectral band (IFitzpatrickl Il999l ) . 



5.1.1 Swift-XRT 



All the Swift-X-ray Telescope (XRT; iBurrows et al.l 120051 ) Windowed Timing observa- 



tions carried out from MJD 54905 to 55044 were used for the analysis: this amounts to a total 
of 41 observations performed within this dedicated multi-instrument effort to study Mrk501. 
The XRT data set was first processed with the XRTDAS software package (v. 2. 5.0) developed 
at the ASI Science Data Center (ASDC) and distributed by HEASARC within the HEASof t 
package (v.6.7). Event files were calibrated and cleaned with standard filtering criteria with 
the xrtpipeline task using the latest calibration files available in the Swift CALDB. The 
individual XRT event files were then merged together using the XSELECT package and the 
average spectrum was extracted from the summed event file. Events for the spectral analysis 
were selected within a circle of 20-pixel (~ 47arcsec) radius centered at the source position 
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and enclosing about 95% of the point-spread function (PSF) of the instrument. The back- 
ground was extracted from a nearby circular region of 40-pixel radius. The source spectrum 
was binned to ensure a minimum of 20 counts per bin to utilize the \ 2 minimization fitting 
technique. The ancillary response files were generated with the xrtmkarf task applying cor- 
rections for the PSF losses and CCD defects using the cumulative exposure map. The latest 
response matrices (v. Oil) available in the Swift CALDB were used. 

The XRT average spectrum in the 0.3 — lOkeV energy band was fitted using the 
XSPEC package. We adopted a log-parabolic model for the photon flux spectral density 
jMassaro et al.ll2004al Jbh of the form \og[T(E)] = logK-a log[£/keV] -b log 2 [£/keV], with 
an absorption hydrogen-equivalent colu mn density fixed to the Galactic value in the direction 
of the source, namely 1.56 x 10 cm" " 2 kalberla et al.ll2005h . This model provided a good 
description of the observed spectrum, with the exception of the 1.4 — 2.3 keV energy band 
where spectral fit residuals were present. These residuals are due to known XRT calibration 
uncertainties (SWIFT-XRT-CALDB-12)E_ and hence we decided to exclude the 1.4-2.3keV 
energy band from the analysis. In addition, we had to apply a small energy offset (~ 40 eV) 
to the observed energy spectrum. The origin of this correction is likely to be CCD charge 
traps generated by radiation and high-energy proton damage (SWIFT-XRT-CALDB-12), 
which affects mostly the lowest energies (first one or two bins) in the spectrum. The result- 
ing spectral fit gave the following parameters: K = (3.41 ± 0.03) x 10 -2 phcm -2 s _1 keV -1 , 
a = 1.96±0.04, and b = 0.308±0.010. The XRT SED data shown in Figure[8]were corrected 
for the Galactic absorption and then binned into 10 energy intervals. 



5.1.5. RXTE-PCA 



The Rossi- X-ray Timing Explorer (RXTE; iBradt et al.l Il993[ ) satellite performed 29 
pointing observations of Mrk501 during the time interval MJD 54905 and 55044. These 
observations amount to a total exposure of 52 ks, which was requested through a dedicated 
Cycle 13 proposal to provide X-ray coverage for our campaign. We did not find a significant 
signal in the RXTE-HEXTE data and hence we only report on the data from RXTE-PCA, 
which is the main pointing instrument on board RXTE. The data analysis was performed 
using FT00LS v6.5 and following the procedures and filtering criteria recommended by the 



RXTE Guest Observer Facility 14 ! after September 2007. In particular, the observations 



13 http://lieas arc. gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/SWIFT-XRT-CALDB-09_vl2.pdf 



14- 



http : //www . universe . nasa . gov/xrays/ programs/rxte/pca/doc/bkg/bkg-2007-saa/ 
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were filtered following the conservative procedures for faint sources^: Earth elevation angle 
greater than 10°, pointing offset less than 0.02°, time since the peak of the last SAA (South 
Atlantic Anomaly) passage greater than 30 minutes, and electron contamination less than 
0.1. For further details on the analysis of faint sources with RXTE, see the online Cook 
BoolJ^I. In the data analysis, in order to increase the quality of the signal, only the first 
xenon layer of PCU2 was used. We used the package pcabackest to model the background 
and the package saextrct to produce spectra for the source and background files and the 
script^ pcarsp to produce the response matrix. 

The PCA average spectrum in the 3 — 28 keV energy band was fitted using the XSPEC 
package with a single power-law function log [J 7 (.£7)] = log A" — a log[E/keV) with a constant 
neutral hydrogen column density N-r fixed at the G alactic value in the direction of the 
source, namely 1.56 x 10 cm" " 2 fcalberla et alil2005h . However, since the PCA bandpass 



starts at 3 keV, the value used for Ah does not significantly affect our results. The resulting 
spectral fit provided a good representation of the data for the following parameters: K = 
(4.34 ± 0.11) x lO^phcm^s^keV" 1 , and a = 2.28 ± 0.02. The PCA average spectrum 
obtained using 23 energy bins is shown in Figure El 



5.1.6. Swift- BAT 



The Swift-Burst Alert Telescope (BAT; lBarthelmy et al. I l2005h analysis results presented 
in this paper were derived with all the available data during the time interval MJD 54905 



and 55044. The spectrum was extracted following the recipes presented in lAjello et al. 



(120081 . l2009bl ). This spectrum is constructed by weighted averaging of the source spectra 
extracted from short exposures (e.g., 300s) and is representative of the averaged source 
emission over the time range spanned by the observations. These spectra are accurate to 
the mCrab level and the reader is referred to lAjello et al.l (j2009at ) for more details. The 
Swift-BAT spectrum is consistent with a power-law function with normalization parameter 
K = 0.24 ± O.iephcm^s^keV" 1 and photon index a = 2.8 ± 0.4. 



15 The average net count rate from Mrk501 was about 7ct/s/pcu (in the energy range 3 — 20keV) with 
flux variations typically much smaller than a factor of two. 

le http : / /heasarc . gsf c . nasa . gov/docs/xte/recipes/ cook_book .html 



17 The CALDB files are located at http://heasarc.gsfc.nasa.gov/FTP/caldb 
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5.1.7. MAGIC 



MAGIC is a system of two 17m-diameter IACTs for very high energy 7-ray astron- 
omy located on the Canary Island of La Palma, at an altitude of 2200 m above sea level. 
At the time of the observation, MAGIC-II, the new second telescope of the current array 
system, was still in its commissioning phase so that Mrk501 wa s observed in stan d-alone 



mode by MAGIC-I, which is in scientific operation since 2004 ([Albert et al.l 120081 ). The 
MAGIC telescope monitored the VHE activity of Mrk501 in the framework of the orga- 
nized multifrequenc y campaign. The observations were performed in the so-called "wobble" 



mode ( ]Daumlll997l ). In order to have a low energy threshold, only observations at zenith 
angles less than 35° were used in this analysis. Bad weather and a shut-down for a sched- 
uled hardware system upgrade during the period MJD 54948-54960 (April 27 - May 13) 
significantly reduced the actual amount of observing time compared to what had initially 



m 



Albert et al. 


(2008 


) and 


Aliu et al. 


(2009) 



to a total of 16.2 hours. The preliminary reconstructed photon fluxes for the individual 
observations gave an average activity of about 30% the flux of the Crab Nebula, with small 
(typically much less than a factor of two) flux variations. The derived spectrum was unfolded 
to correct for the eff ects of the limited energy resolution of the detector and of possible bias 
( 1 Albert et al.ll2007bl ). The resulting spectrum was fitted satisfactorily with a single power- 
law function of the form log[J-"(£/)] = logK — a log[£J/TeV], giving normalization parameter 
K = (0.90 ± 0.05) x 10- n phcm~ 2 s' 1 TeV' 1 and photon index a = 2.51 ± 0.05. 



5.1.8. VERITAS 



VERITAS is a state-of-the-art TeV 7-ray observatory consisting of four 12 m-diameter 
IACTs. VERITAS is located at the basecamp of the F.L. Whipple Observatory in southern 
Arizona, USA, at an alti tude of 1250 m aboy e sea level, and the system has been fully op- 
erational since fall 2007 ( lAcciari et ajj |2010aj ). VERITAS observed Mrk501 as part of the 
long-term monitoring campai gn between March and June of 2009. The observations were 
performed in "wobble" mode (jDaumll 19971 ) at relative ly low zenith angle ( < 40°). These data 
were analyzed following the prescription reported in lAcciari et al.l ( 120081 ). After removal of 
data runs with poor observing conditions, a total of 9.7 hours of good quality data was 
obtained between MJD 54907 and MJD 55004. Due to the long-term nature of these obser- 
vations, several factors had to be taken into account when analyzing the data. The initial 
portion of the campaign includes data taken under standard 4-telescope operating condi- 
tions. Two nights of data were taken with only two operational telescopes due to technical 
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difficulties. For the latter portion of the campaign, data were taken over several nights with 
three operational tele scopes because one of the tel escopes was being relocated as part of an 



upgrade to the array (IPerkins &: Maier et al.ll2009l ). The effective collection areas for the ar- 



ray in these three configurations were calculated using Monte Carlo simulations of extensive 
air showers passed through the analysis chain with detector configurations corresponding to 
the respective data-taking conditions. 

An initial analysis of the VHE activity showed an increase in the flux by a factor of about 
five during MJD 54953-54956. Because of the large difference in the VHE flux, we decided 
to analyze this 3-day data set (corresponding to a "flaring" state of Mrk501) separately from 
the rest of the collected data ( "non-flaring" ) . The "flaring" epoch consists of 2.4 h of data 
taken during MJD 54953-54956. The "non-flaring" epoch consists of 7.3 h of data taken 
during the remaining portion of the campaign. The spectra from these two data sets were 
each fitted with a single power-law function of the form log[J r (-E')] = logK — a log[E/TeV}. 
The resulting fit parameter values are K = (4.17 ± 0.24) x 10~ n phcm -2 s _1 TeV -1 with 
a = 2.26 ±0.06 for the "flaring" state, and K = (0.88 ± 0.06) x 10~ n ph cm" 2 s" 1 TeV" 1 with 
photon index a = 2.48 ± 0.07 for the "non-flaring" state. 



Table 1. List of instruments participating in the multifrequency campaign and used in the construction of the SED 

in FigureE] 



Instrument / observatory 



Energy range covered 



Web page 



MAGIC 

VERITAS 

Whipple" 

Fermi-LAT 

Swift-BAT 

RXTE-PCA 

Swift-XRT 

Sra/i-UVOT 

Abastumani (through GASP-WEBT program) 
Lulin (through GASP-WEBT program) 

Roque de los Muchachos (KVA) (through GASP-WEBT program) 

St. Petersburg (through GASP-WEBT program) 

Talmassons (through GASP-WEBT program) 

Valle d'Aosta (through GASP-WEBT program) 

GRT 

MitSume 

ROVOR 

Campo Imperatore (through GASP-WEBT program) 

OAGH 

WIRO 

SMA 

VLBA 

Noto 

Metsahovi (through GASP-WEBT program) 
VLBA (through MOJAVE program) 
OVRO 
Medicina 

UMRAO (through GASP-WEBT program) 
RATAN-600 

EfTclsberg (through F-GAMMA program) 



0.12-5.8 TeV 
0.20-5.0 TeV 
0.4-1.5 TeV 
0.1-400 GeV 
14-195 keV 
3-28 keV 
0.3-9.6 keV 

V, B, U, UVW1, UVM2, UVW2 

R band 

R band 

R band 

R band 

R band 

R band 

V, R, B bands 

g, Re, Ic bands 

B, R, V, I bands 

H, J, K bands 

H, J, K bands 

J, K bands 

225 GHz 

4.8, 8.3, 15.4, 23.8, 43.2 GHz 

8.4, 43 GHz 

37 GHz 

15 GHz 

15 GHz 

8.4, 22.3 GHz 

4.8, 8.0, 14.5 GHz 

2.3, 4.8, 7.7, 11.1, 22.2 GHz 

2.6, 4.6, 7.8, 10.3, 13.6, 21.7, 31 GHz 



http 
http 
http 
http 




http 
http 
http 
http 
http 



http 



http 



http 



http 
http 
http 



http 
http 
http 
http 
http 
http 



http 



http 



//wwwmagic .mppmu.mpg.de/ 
//veritas . sao . arizona.edu/ 

//veritas . sao . arizona. edu/content/blogsect ion/6/40/ 
//www-glast . Stanford. edu/ index .html 



//heasarc . gsfc .nasa.gov/docs/swif t/swiftsc .html/ 
//heasarc . gsfc .nasa.gov/docs/xte/rxte .html 
//heasarc . gsfc .nasa.gov/docs/swif t/swiftsc .html 
//heasarc . gsfc .nasa.gov/docs/swif t/swiftsc .html 
//www. oato . inaf . it/blazars/webt/ 



//www. oato . inaf . it/blazars/webt/ 



//www. oato . inaf . it/blazars/webt/ 

//www. oato . inaf . it/blazars/webt/ 

//www. oato . inaf . it/blazars/webt/ i 

//www. oato . inaf . it/blazars/webt/ 

//asd. gsf c .nasa. gov/Takanori . Sakamoto/GRT/index^html 
//www.hp.phys . titech. ac . jp/mitsume/index .html 
//rovor . byu . edu/ 



//www. oato . inaf . it/blazars/webt/ 



//astro . inaoep.mx/en/observatories/oagh/ 
//physics .uwyo . edu/- chip/wiro/wiro .html 
//smal . sma . hawaii . edu/ 
//www . vlba . nrao . edu/ 
//www. noto . ira. inaf . it/ 
//www. metsahovi . f i/ 
//www. physics .purdue . edu/MOJAVE/ 
//www. astro . caltech . edu/ovroblazars 
//www.med. ira. inaf . it/index_EN.htm 
//www. oato . inaf . it/blazars/webt/ 



//wO . sao . ru/ratan/ 



//www . mpif r-bonn . mpg . de/div/ef f elsberg/index_e . html/ 
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Fig. 1.— Left: Fermz-LAT 7-ray flux in the energy range 0.3 — 400 GeV (top panel) 
and spectral photon index from a power-law fit (bottom panel) for Mrk501 for 30-day time 
intervals from 2008 August 5 (MJD 54683) to 2009 November 27 (MJD 55162). Vertical bars 
denote la uncertainties and the horizontal bars denote the width of the time interval. The 
red dashed line and the red legend show the results from a constant fit to the time interval 
MJD 54862-54982, while the black dashed line and black legend show the results from a 
constant fit to the entire 480-day data set. Right: Scatter plot of the photon index vs flux 
values. 
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Fig. 2. — Multifrequency light curves of Mrk501 with 30-day time bins obtained with 3 
all-sky- monitoring instruments: RXTE-ASM (2 — lOkeV, first from the top); Swift-B AT 
(15 — 50keV, second) and Fermi-LAT for two different energy ranges (0.2 — 2GeV, third, 
and > 2 GeV, fourth). The light curves cover the period from2008 August 5 (MJD 54683) to 
2009 November 27 (MJD 55162). Vertical bars denote la uncertainties and horizontal bars 
show the width of the time interval. The horizontal dashed lines and the legends (for all the 
plots) show the results from a constant fit to the entire 480-day data set. 
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Fig. 3. — Fractional variability parameter for 16 months data (2008 August 5 - - 2009 
November 27) from 3 all-sky- monitoring instruments: RXTE- ASM (2 — lOkeV); Swift- 
BAT (15 - 50keV) and Fermi-LAT (two energ y ranges 0.2 - 2 GeV a nd 2 - 300 GeV). The 
fractional variability was computed according to lVaughan et al.l (120031 ) using the light curves 
from Figure[2l Vertical bars denote la uncertainties and horizontal bars indicate the width of 
each energy bin. The horizontal dashed line and the legend show the results from a constant 
fit. 
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Fig. 4. — Spectral energy distribution for Mrk501 from Fermi-LAT during the period from 
2008 August 5 (MJD 54683) to 2009 November 27 (MJD 55162). The black line depicts the 
result of the unbinned likelihood power-law fit, the red contour is the 68% uncertainty of 
the fit, and the black data points show the energy fluxes in differential energy ranges. The 
legend reports the results from the unbinned likelihood power-law fit in the energy range 
0.3 - 400 GeV. 
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MJD 54892-54922; TS = 171 
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Fig. 5. — Spectral energy distribution for Mrk501 from Fermi-LAT for six 30-day time 
intervals: MJD 54832-54862 (top left), MJD 54862-54892 (top right), MJD 54892-54922 
(middle left), MJD 54922-54952 (middle right), MJD 54952-54982 (bottom left) and MJD 
54982-55012 (bottom right). In all the panels, the black line depicts the result of the 
unbinned likelihood power-law fit, the red contour denotes the 68% uncertainty of the fit, 
and the black data points show the energy fluxes computed for differential energy ranges. 
The blue arrows denote 95% upper limits, which were computed for the differential energy 
ranges with a signal of TS < 4 or less than two photons predicted by the analysis model for 
Mrk501. The legend reports the results from the unbinned likelihood power- law fit in the 
energy range 0.3 — 400 GeV. 
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Fig. 6. — Time and energy coverage during the multifrequency campaign. For the sake of 
clarity, the minimum observing time displayed in the plot was set to half a day. 



Note. — The energy range shown in column two is the actual energy range covered during the Mrk501 observations, and not the instrument's nominal energy range, which 
might only be achievable for bright sources and excellent observing conditions. 



Note. — (a) The Whipple spectra were not included in Figure[8] See text for further comments. 



Co 
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5.2. Fermi-hAT Spectra During the Campaign 

The Mrk501 spectrum measured by Fermi-LAT, integrated during the time interval 
of the multifrequency campaign, is shown in the panel (b) of Figured The spectrum can 
be described by a power-law function with photon index 1.74 ± 0.05. The flux data points 
resulting from the analysis in differential energy ranges are within 1 — 2 a of the power- 
law fit result; this is an independent indication that a single power-law function is a good 
representation of the spectrum during the multifrequency campaign. On the other hand, 
the shape of the spectrum depicted by the differential energy flux data points suggests the 
possibility of a concave spectrum. As it was discussed in §3] and §H (see Figures[T] and 
|5J), Mrk501 showed substantial spectral variability during the time period covered by the 
multifrequency campaign, with some 30-day time intervals characterized by relatively soft 
spectra (photon index ~ 2 for the 30-day intervals MJD 54892-54922 and MJD 54922-54952) 
and others by relatively hard spectra (photon index ~ 1.6 for the 30-day intervals MJD 
54952-54982, MJD 54982-55012 and MJD 55012-55042). The panel (b of Figured presents 
the average spectrum over those time intervals, and hence it would not be surprising to see 
two slopes (instead of one) in the spectrum. In order to evaluate this possibility, a broken 
power-law fit was applied, yielding indices of 1.86 ± 0.08 and 1.44 ± 0.14 below and above a 
break energy of 10 ± 3 GeV, respectively. The likelihood ratio of the broken power law and 
the power law is 2.2. Given that the broken power law has two additional degrees of freedom, 
this indicates that the broken power law is not statistically preferred over the single power 
law function. 

For comparison purposes we also computed the spectra for time intervals before and 
after the multifrequency campaign (MJD 54683-54901 and MJD 55044-551620 These two 
spectra, shown in the panel (a) and (ac of Figured can both be described satisfactorily by 
single power-law functions with photon indices 1.82 ± 0.06 and 1.80 ± 0.08. Note that the 
two spectra are perfectly compatible with each other, which is consistent with the relatively 
small flux/spectral variability shown in Figures!]] and |2] for those time periods. 

5.3. The Average Broadband SED During the Campaign 

The average broadband SED of Mrk501 resulting from our 4.5-month-long multifre- 
quency campaign is shown in FigureEJ The TeV data from MAGIC and VERITAS have been 



18 Tcchnical problems prevented the scientific operation of the Fermi-LAT instrument during the interval 
MJD 54901-54905. 
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correc ted for the absorption in the EBL using the particular EBL model bvlFranceschini et al. 



(200 81). The correc t ions given by the other low-EBL-level models (IKneiske et al 



2004 



Gilmore et al.l 120091 ; iFinke et al.l 120101 ) are very similar for the low redshift of Mrk501 
(z = 0.034). The attenuation factor at a photon energy of 6 TeV (the highest energy detected 
from Mrk501 during this campaign) is in the range e~ Tl "< ~ 0.4 — 0.5, and smaller at lower 
energies. 

During the campaign, as already noted above, the source did not show large flux vari- 
ations like those recorded by EGRET in 1996, or those measured by X-ray and TeV in- 
struments in 1997. Nevertheless, significant flux and spectral variations at 7-ray energies 
occurred in the time interval MJD 54905-55044. The largest flux variation during the cam- 
paign was observed at TeV energies during the time interval MJD 54952.9-54955.9, when 
VERITAS measured a flux about five times higher than the average one during the cam- 
paign. Because of the remarkable difference with respect to the rest of the analyzed exposure, 
these observations were excluded from the data set used to compute the average VERITAS 
spectrum for the campaign; the three-day "flaring-state" spectrum (2.4 hours of observation) 
is presented separately in FigureO Such a remarkable flux enhancement was not observed 
in the other energy ranges and hence Figure[S] shows only the averaged spectra for the other 
instrument^]. 

The top panel in Figured shows a zoom of the high-energy bump depicted in FigureEJ 
The last two energy bins from Fermi (60 — 160 and 160 — 400 GeV) are systematically above 
( l-2cr ) the measured/extrapolated spectrum from MAGIC and VERITAS. Even though this 
mismatch is not statistically significant, we believe that the spectral variability observed 
during the 4.5 month long campaign (see £]Hand §5.21) could be the origin of such a difference. 
Because Fermi-LAT operates in a survey mode, Mrk501 is constantly monitored at GeV 
energies^!, while this is not the case for the other instruments which typically sampled the 
source during <1 hour every 5 days approximately. Moreover, because of bad weather or 
moonlight conditions, the monitoring at the TeV energies with Cherenkov telescopes was even 
less regular than that at lower frequencies. Therefore, Fermi-LAT may have measured high 
activity that was missed by the other instruments. Indeed, the 2.4-hour high-flux spectrum 
from VERITAS depicted in Figure[S] (which was obtained during the 3-day interval MJD 
54952.9-54955.9) demonstrates that, during the multifrequency campaign, there were time 
periods with substantially (factor of five) higher TeV activity. It is possible that the highest- 



19 The MAGIC telescope did not operate during the time interval MJD 54948-54965 due to a drive system 
upgrade. 

20 During every three hours of Fermi operation, Mrk501 is in the LAT field of view for about 0.5 hour. 
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energy LAT observations (>50 GeV) include high TeV flux states which occurred while the 
IACTs were not observing. 

If the flaring activity occurred only at the highest photon energies, then the computed 
Fermi-LAT flux (>0.3 GeV) would not change very much and the effect might only be visible 
in the measured power-law photon index. This seems to be the case in the presented data set. 
As was shown in Figure^ the 30-day intervals MJD 54922-54952 and MJD 54952-54982 have 
photon fluxes above 0.3 GeV of (3.9±0.6) x 10~ 8 ph cm~ 2 s" 1 and (3.6±0.5) x 10~ 8 ph cm^s" 1 , 
while their photon indices are 2.10 ± 0.13 and 1.63 ± 0.09, respectively. Therefore, the 
spectral information (together with the enhanced photon flux) indicates the presence of 
flaring activity at the highest 7-ray energies during the second 30-day time period. Besides 
the factor ~ 5 VHE flux enhancement recorded by VERITAS and Whipple at the beginning 
of the time interval MJD 54952-54982, MAGIC and Whipple also recorded a factor ~ 2 VHE 
flux enhancement at the end of t his 30-day time interval (see preliminary fluxes reported in 



Panequdl2009i ; iPichel et al.ll2009i ). This flux enhancement was measured for the time interval 
MJD 54975-54977, but there were no VHE measurements during the period MJD 54970.5- 
54975.0. Thus, the average Fermi-LAT spectrum could have been affected by elevated VHE 
activity during the 30-day time interval MJD 54952-54982, which was only partly covered 
by the IACTs participating in the campaign. 

For illustrative purposes, in the bottom panel of Figured we show separately the Fermi- 
LAT spectra for the 30-day time interval MJD 54952-54982 (high photon flux and hard 
spectrum), and for the rest of the campaign. It is interesting to note that the Fermi-LAT 
spectrum without the 30-day time interval MJD 54952-54982 (blue data points in the bottom 
panel of Figure[9]) agrees perfectly with the VHE spectrum measured by IACTs. We also want 
to point out that the power-law fit to the Fermi-LAT spectrum without the 30-day interval 
MJD 54952-54982 gave a photon flux above 0.3 GeV of (2.62 ± 0.25) x lO^phcm" 2 s" 1 
with a photon index of 1.78 ± 0.07, which is statistically compatible with the results for the 
power-law fit to the Fermi-LAT data from the entire campaign (see panel (b) in Figure[7]). 
As discussed above, the flaring activity occurred mostly at the highest energies, where the 
(relatively) low photon count has little impact on the overall power-law fit performed above 
0.3 GeV. 

This is the most complete quasi-si multaneous SEP ev er collected for Mrk501, or for 



any other TeV-emitting BL Lac (see also lAbdo et al.ll2010dl . in preparation). At the highest 



energies, the combination of Fermi and MAGIC/VERITAS data allows us to measure, for 
the first time, the high-energy bump without any spectral gap. The low-energy spectral 
component is also very well characterized with Swift-XJVOT , Swift-XKT and RXTE-PCA 
data, covering the peak of the synchrotron continuum. The only (large) region of the SED 
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with no available data corresponds to the photon energy range 200 keV — 100 MeV, where 
the sensitivity of current instruments is not good enough to detect Mrk501. It is worth 
stressing that the excellent agreement in the overlapping energies among the various instru- 
ments (which had somewhat different time coverage) indicates that the collected data are 
representative of the truly average SED during the multi-instrument campaign. 



6. Modeling the Spectral Energy Distribution of Mrk501 



The simultaneous broadband data set resulting from the multifrequency campaign re- 
ported above offers an unprecedented opportunity to model the emission of an archetypal 
TeV blazar in a more robust way than in the past. It is widely believed that the radio-to-7- 
ray emission of the BL Lac class of AGN is produced predominantly via the synchrotron and 
synchrotron self-Compton (SSC) processes, and hence the homogeneous one-zone approx- 
imation of the SSC scenario is the simplest model to consider. Here we therefore adopt 
the 'standard' one-zone SSC model, which has had moderate success in accounting for 
the spectral and temporal properties of the TeV-emitting BL Lacs analyzed so far (e.g., 



Finke et al.l 120081 ; iGhisellini et al.l l2009al and references therein). We also note that one- 
zone SSC analyses have been widely applied before to the particular case of Mrk501 (e.g. 



Bednarek &: 



Albert et al 



rotheroe 1999; Katarzvhski et al.ll200ll ; lTavecchio et al.ll200ll ; iKino et al.ll2002 



2007af ). However, it is important to stress that the modeling results from the 



previous works related almost exclusive ly to the high-activity state of Mrk501. In the more 
recent work by lAnderhub et al.l (120091 ) the source was studied also during its low-activity 
state, yet the simultaneous observations used in the modeling covered only the X-ray and 
TeV photon energies. In this paper we study Mrk501 during a relatively low activity state, 
and the modeling is applied to a more complete broadband SED extending from radio to 
TeV energies, including the previously unavailable GeV data from Fermi. This constitutes 
a substantial difference with respect to previous works. The resulting constraints on the 
physical parameters of the source, together with several limitations of the applied scenario, 
are discussed further down in the next sections. 

We want to note that modeling of the average blazar SED based on a scenario assum- 
ing steady-state homogenous emission zone could be an over-simplification of the problem. 
The blazar emission may be produced in an inhomogeneous region, involving stratification 
of the emitting plasma both along and accross a relativistic outflow. In such a case, the 
observed radiative output of a blazar could be due to a complex superposition of different 
emission zones characterized by very different parameters and emission properties. Some 
first attempts to approach this problem in a more quantitative way have been already dis- 
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cussed in the literatu re (e.g. 



Ghisellini et al. 


2005; 


Katarzvnski et al. 


2008; 


Graff et al. 


2008 



Giannios et al.l 120091 ) . The main drawback of the proposed models, however, is the increased 



number of free parameters (over the simplest homogeneous one- zone scenario), what reduces 
considerably the predictive power of the modeling. That is particularly problematic if a 
"limited" (in a time and energy coverage) dataset is considered in the modeling. Only a 
truly simultaneous multifrequency dataset covering a large fraction of the available elec- 
tromagnetic spectrum and a wide range of timescales — like the one collected during this 
and future campaigns which will be further exploited in forthcoming publications — will 
enable to test such more sophisticated and possibly more realistic blazar emission models in 
a time-deoendent manner. 



6.1. SSC Modeling 

Let us assume that the emitting region is a homogeneous and roughly spherically sym- 
metric moving blob, with radius R and comoving volume V ~ (4tt/3)R 3 . For this, we 
evaluate the comoving synchrotron and synchrotron self-Compton emissivities, u'j'^,, assum- 
ing isotropic distributions of ultrarelativistic electrons and synchrotron photons in the rest 
frame of the emitting region. Thus, we use the exact synchrotron and inverse- Compton 
kernels (the latter on e valid in b oth Thomson and Kl e in-Nis hina regimes), as given in 



Crusius fc Schlickeiserl ( 119861 ) and iBlumenthal fc Gouldl (Il970l ). respectively. The intrin- 



sic monochromatic synchrotron and SSC luminosities are then v'L' v , = An V h>'j' u , , while the 
observed monochromatic flux densities (measured in erg cm -2 s -1 ) can be found as 

vF v = W L 'v<]v>=v(l+z)/5 - [ V %']v'=u(l+z)/5 , (2) 

where 5 is the jet Doppler factor, z = 0.034 is the source redshift, and = 142 Mpc is 
the luminosity distance to Mrk501. In order to evaluate the comoving emissivities u'j' v ,, the 
electron energy distribution ti^Ct) nas ^° be specified. For this, we assume a general power- 
law form between the minimum and maximum electron energies, 7 mm and 7 ma x, allowing 
for multiple spectral breaks in between, as well as for an exponential cut-off above 7 ma x- In 
fact, the broadband data set for Mrk501 requires two different electron break energies, and 
hence we take the electron energy distribution in a form 

{7~ S1 for 7mm < 7 < Ibr, 1 

1~ S2 for 7br,i < 7 < 7&r,2 , (3) 

7~ S3 exp [-7/7 mox ] for lbr,2 < 7 
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with the normalization expressed in terms of the equipartition parameter (the ratio of the 
comoving electron and magnetic field energy densities), namely 

= I im e c 2 n' e {-i) d-f 

The measured SED is hardly compatible with a simpler form of the electron distribution 
with only one break and an exponential cutoff. However, some smoothly curved spec- 
tral shape might perhaps be an alternative repr esentation of the electron spectrum (e.g., 



Stawarz &: Petrosianl 120081 ; iTramacere et al.ll2Q09l ) 



/ max i 



The model adopted is thus characterized by four main free parameters (B, R, 5, and r\ e 
plus seven additional ones related to the electron energy distribution (7 m j„, 7& r ,i, 76r,2> 7r? 
si, S2, and S3). These seven additional parameters are determined by the spectral shape of the 
non-thermal emission continuum probed by the observations, predominantly by the spectral 
shape of the synchrotron bump (rather than the inverse-Compton bump), and depend only 
slightly on the particular choice of the magnetic field B and the Doppler factor 5 within the 
allowed ranged . There is a substantial degeneracy regarding the four main free parameters: 
the average emission spectrum of Mrk501 may be fitted by different combinations of B, R, 
5, and rj e with little variation in the shape of the electron energy distribution. Note that, 
for example, [vF v ] syn oc R 3 r) e , but at the same time [zaFj,] ssc oc R 4 r] 2 . We can attempt to 
reduce this degeneracy by assuming that the observed main variability timescale is related 
to the size of the emission region and its Doppler factor according to the formula 

(l + z)R 

(5) 



CO 

The multifrequency data collected during the 4.5-month campaign (see §5]) allows us to study 
the variability of Mrk501 on timescales from months to a few days. We found that, during 
this time period, the multifrequency activity varied typically on a timescale of 5 — 10 days, 
with the exception of a few particular epochs when the source became very active in VHE 
7-rays, and flux variations with timescales of a day or shorter were found at TeV energies. 
Nevertheless, it is important to stress that several authors concluded in the past that the 
domina nt emission site of M rk501 is characterized by variability timescales longer than one 



day (see lKataoka et al.ll200lL for a comprehensive study of the Mrk501 variability in X-rays), 
and that the power in the intraday flickering of this source is small, in agreement with the 
results of our campaign. Nevertheless, one should keep in mind that this object is known 
for showing sporadic but extreme changes in its activity that can give flux variations on 



21 For example, for a given critical (break) synchrotron frequency in the observed SED, the corresponding 
electron break Lorentz factor scales as 7t, r oc 1/ \/B S. 
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timescales as short as a few minutes ([Albert et all l2007al ). In this work we aim to model 



the average/typical behaviour of Mrk501 (corresponding to the 4.5-month campaign) rather 
than specific/short periods with outstanding activity, and hence we constrained the minimum 
(typical) variability timescale t var in the model to the range 1 — 5 days. 

Even with t var fixed as discussed above, the reconstructed SED of Mrk501 may be fitted 
by different combinations of B, R, 5, and r] e . Such a degeneracy bet ween the main mode l 



parameters is an inevitable feature of the SSC modeling of blazars (e.g-. lKataoka et al.lll999l ). 
and it is therefore necessary to impose additional constraints on the physical parameters of 
the dominant emission zone. Here we argue that such constraints follow from the requirement 
for the electron energy distribution to be in agreement with the one resulting from the 
simplest prescription of the energy evolution of the radiating electrons within the emission 
region, as discussed below. 

The idea of separating the sites for the particle acceleration and emission processes is 
commonly invoked in modeling different astrophysical sources of high-energy radiation, and 
blazar jets in particular. Such a procedure is not always justified, because interactions of 
ultrarelativistic particles with the magnetic field (leading to particle diffusion and convec- 
tion in momentum space) are generally accompanied by particle radiative losses (and vice 
versa). On the other hand, if the characteristic timescale for energy gains is much shorter 
than the timescales for radiative cooling (t' rad ) or escape (t' esc ) from the system, the particle 
acceleration processes may be indeed approximated as being 'instantaneous', and may be 
modeled by a single injection term Q("f) in the simplified version of the kinetic equation 



dt ^7 



7^(7) 



ne(7) +Q(7) (6) 



t 



describing a very particular scenario for the energy evolution of the radiating ultrarelativistic 
electrons. 

It is widely believed that the above equation is a good approximation for the energy evo- 
lution of particles undergoing diffusive ('first-order Fermi') shock acceleration, and cooling 
radiatively in the downstream region of the shock. In such a case, the term speci- 
fies the energy spectrum and the injection rate of the electrons freshly accelerated at the 
shock front and not affected by radiative losses, while the escape term corresponds to the 
energy-independent dynamical timescale for the advection of the radiating particles from the 
downstream region of a given size R, namely t' esc ~ t' d ~ R/c. The steady-state electron 
energy distribution is then very roughly n'J^y) ~ t' d Qij) below the critical energy for which 
Kadi"?) = ^'dyni an d n ' e (l) ~ Kadil) Q(l) above this energy. Note that in the case of a power- 
law injection Q(j) oc 7 _s and a homogeneous emission region with dominant radiative losses 
of the synchrotron type, t' rad (j) oc 7 -1 , the injected electron spectrum is expected to steepen 
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by As = 1 above the critical 'cooling break' energy. This provides us with the additional 
constraint on the free model parameters for Mrk501: namely, we require that the position 
of the second break in the electron energy distribution needed to fit the reconstructed SED, 
7b r 2, should correspond to the location of the cooling break for a given chosen set of the 
model free parameters. 

Figure[10] (black curves) shows the resulting SSC model fit (summarized in Table 2) to 
the averaged broadband emission spectrum of Mrk501, which was obtained for the following 
parameters: B = 0.015 G, R = 1.3 X 10 17 cm, 5 = 12, rj e = 56, ^ min = 600, 7^1 = 4 x 10 4 , 
Ibr, 2 = 9 x 10 5 , jmax = 1-5 x 10 7 , Si = 2.2, s 2 = 2.7, and S3 = 3.65. The overall good 
agreement of the model with the data is further discussed in §6.21 Here, we note that, for 
these model parameters, synchrotron self- absorption effects are important only below 1 GHz, 
where we do not have observations^]. We also emphasize that with all the aforementioned 
constraints and for a given spectral shape of the synchrotron continuum (including all the 
data points aimed to be fitted by the model, as discussed below), and thus for a fixed spectral 
shape of the electron energy distribution (modulo critical electron Lorentz factors scaling as 
oc 1/VB 5), the allowed range for the free parameters of the model is relatively narrow. 
Namely, for the variability timescale between 1 and 5 days, the main model parameters may 
change within the ranges R ~ (0.35 - 1.45) x 10 17 cm, 5 ~ 11 — 14, and B ~ 0.01 - 0.04 G. 
The parameter 7] e depends predominantly on the minimum Lorentz factor of the radiating 
electrons. Hence, it is determined uniquely as r] e ~ 50 with the sub-mm flux included in 
the fitted dataset. Only with a different prescription for the spectral shape of the electron 
energy distribution could the main free parameters of the model be significantly different 
from those given above. 

Despite the absence of any fast variability during this multifrequency campaign (apart 
from the already discussed isolated 3-day-long flare), Mrk501 is known for the e xtremely 



rapid flux changes at the highest observed photon energies (e.g. lAlbert et al.ll2007af ). Hence, 
it is interesting to check whether any shorter than few-day-long variability timescales can be 
accommodated in the framework of the simplest SSC model applied here for the collected 
dataset. In order to do that, we decreased the minimum variability time scale by one order of 
magnitude (from 4 days to 0.4 days), and tried to model the data. A satisfactory fit could be 
obtained with those modified parameters, but only when we relaxed the requirement for the 
electron energy distribution to be in agreement with the one following from the steady-state 
solution to Equation [6j and in particular the resulting constraint for the second break in 



The turnover frequency related t o the synchrotron self-absorption may be evaluated using the formulae 

60 MHz, 



given in iGhisellini fc Svenssonl ([19911 ) and the parameter values from our SSC model fit as v' ss 
which in the observer frame reads v ssa — 8v' ssa /(l + z) ~ 0.7 GHz. 
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the electron spectrum to be equal the cooling break. This "alternative" model fit is shown 
in Figure[10] (red curves) together with the "best" model fit discussed above. The resulting 
model parameters for the "alternative" fit are B = 0.03 G, R = 2 x 10 16 cm, 5 = 22, r) e = 130, 
lmin = 300, lbr>1 = 3 x 10 4 , 7fcr , 2 = 5 x 10 5 , lmax = 3 x 10 6 , s x = 2.2, s 2 = 2.7, and s 3 = 3.5. 
This particular parameter set — which should be considered as an illustrative one, only 
- would be therefore consistent with a minimum variability timescale of 0.36 days, but at 
the price of much larger departures from the energy equipartition (j] e > 100). The other 
source parameters, on the other hand, would change only slightly (see Table 2). Because 
of the mismatch (by factor ~3) between the location of the cooling break and the second 
break in the electron distribution, we consider this "alternative" fit less consistent with the 
hypothesis of steady-state homogeneous one-zone SSC scenario, which is the framework we 
chose to model the broad-band SED of Mrk501 emerging from the campaign. 



6.2. Notes on the Spectral Data Points 

The low-frequency radio observations performed with single-dish instruments have a 
relatively large contamination from non-blazar emission due to the underlying extended jet 
component, and hence they only provide upper limits for the radio flux density of the blazar 
emission zone. On the other hand, the flux measurements by the interferometric instruments 
(such as VLB A), especially the ones corresponding to the core, provide us with the radio 
flux density from a region that is not much larger than the blazar emission region. 

The radio flux densities from interferometric observations (from the VLBA core) are 
expected to be close upper limits to the radio continuum of the blazar emission component. 
The estimated size of the partially- resolved VLBA core of Mrk 501 at 15 GHz and 43 GHz is ~ 
0.14-0.18mas~ 2.9-3.7 xl0 17 cm (with the appropriate conversion scale 0.67pc/mas). The 
VLBA size estimation is the FWHM of a Gaussian representing the brightness distribution of 
the bl ob, which could be approximated as 0.9 times the radius of a corresponding spherical 
blob ( Marscherl 1983 ). That implies that the size of the VLBA core is only a factor 2-3 
larger than the emission region in our SSC model fit (R = 1.3 x 10 17 cm). Therefore, it is 
reasonable to assume that the radio flux density from the VLBA core is indeed dominated 
by the radio flux density of the blazar emission. Forthcoming multi-band correlation studies 
(in particular VLBA and SMA radio with the 7-rays from Fermi-LAT) will shed light on this 
particular subject. Interestingly, the magnetic field claimed for the partially-resolved radio 
core of Mrk 501 (which has a size of < 0.2 mas) and its sub-mas jet, namely B ~ (10 — 30) mG 



( IGiroletti et al.l l2004j . 120081 ). is in very good agreement with the value emerging from our 



model fits (15 mG), assuring self-consistency of the approach adopted. 
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In addition to this, in the modeling we also aimed at matching the sub-millimeter flux 
of Mrk501, given at the observed frequency of 225 GHz, assuming that it represents the low- 
frequency tail of the optically-thin synchrotron blazar component. One should emphasize in 
this context that it is not clear if the blazar emission zone is in general located deep within 
the millimeter photosphere, or not. However, the broadband variability of luminous blazars 
of the FSRQ type indicates that there is a significant overlap of the blazar zone with a region 



wher e the jet becomes optically thin at millimeter wavelengths (as discussed by lSikora et al. 



20081 . for the particular case of the blazar 3C 454.3). We have assumed that the same holds 



for BL Lac objects. 

The IR/optical flux measurements in the range ~ (1—10) x 10 14 Hz represent the starlight 
of the host galaxy and hence they should be excluded when fitting the non-thermal emission 
of Mrk501. We modelled these data points with the template spectrum of an elliptical galaxy 
instead (inc luding only the dominant stellar component due to the evolved red giants, as 



discussed in lSilva et al.lll998[ ). obtaining a very good match (see the dotted line in Figure lTO]) 
for the bolometric starlight luminosity L sta r — 3 x 10 44 ergs _1 . Such a luminosity is in fact 
expected for the elliptical host of a BL Lac object. The model spectrum of the galaxy falls 
off very rapidly above 5 x 10 14 Hz, while the three UV data points (above 10 15 Hz) indicate a 
prominent, flat power-law UV excess over the starlight emission. Therefore, it is reasonable 
to assume that the observed UV fluxes correspond to the synchrotron (blazar-type) emission 
of Mrk501 and, consequently, we used them in our model fit. However, many elliptical 
galaxies do reveal in their spectra the so-called 'UV upturn', or 'UV excess', whose origin 
is not well known, but which is presumably related to the starlight continuum (most likely 
due to young stars from the residual star-forming activity within the cent ral region of a 



galaxy) rather than to non-thermal (jet-related) emission processes (see, e.g.. lCode fc Welch 



19791 ; lAtlee et al.ll2009l ). Hence, it is possible that the UV data points provided here include 
some additional contamination from the stellar emission, and as such might be considered 
as upper limits for the synchrotron radiation of the Mrk501 jet. 

The observed X-ray spectrum of Mrk 501 agrees very well with the SSC model fit, except 
for a small but statistically significant discrepancy between the model curve and the first two 
data points provided by Swift-XKT, which correspond to the energy range 0.3 — 0.6 keV. As 
pointed out in §5.1} the Swift-XKT data had to be corrected for a residual energy offset which 
affects the lowest energies. The correction for this effect could introduce some systematic 
differences with respect to the actual fluxes detected at those energies. These low-energy 
X-ray data points might be also influenced by intrinsic absorption of the X-ray photons 
within the gaseous envi ronment of Mrk 501 nucleus, as suggested by the earlier studies with 



the ASCA satellite (see lKataoka et al.lll999l ). As a result, the small discrepancy between the 



data and the model curve within the range 0.3 — 0.6 keV can be ignored in the modelling. 
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The agreement between the applied SSC model and the 7-ray data is also very good. In 
particular, the model returns the 7-ray photon index 1.78 in the energy range 0.3-30 GeV, 
which can be compared with the one resulting from the power-law fit to the Fermi-LAT 
data above 0.3 GeV, namely 1.74 ± 0.05. However, the last two energy bins from Fermi 
(60 — 160 and 160 — 400 GeV) are systematically above (2a) the model curves, as well as 
above the averaged spectrum reported by MAGIC and VERITAS. A possible reason for 
mismatch between the average Fermi-LAT spectrum and the one from MAGIC/VERITAS 
was discussed in §5.31 



7. Discussion 

In this section we discuss some of the implications of the model results presented above. 
After a brief analysis of the global parameters of the source resulting from the SSC fits 
( §7.ip . the discussion focuses on two topics. Firstly ( §7.21) . we show that the characteristics 
of the electron energy distribution emerging from our modeling can be used to constrain the 
physical processes responsible for the particle acceleration in Mrk501, processes which may 
also be at work in other BL Lac type objects. Secondly ( §7.3p . we examine the broadband 
variability of Mrk501 in the framework of the model. 



7.1. Main Characteristics of the Blazar Emission Zone in Mrk501 

The values for the emission region size R = 1.3 x 10 17 cm and the jet Doppler factor 
5 = 12 emerging from our SSC model fit give a minimum (typical) variability timescale of 
t var ~ (1 + z)R/cS ~ 4 days, which is consistent with the variabilit y observed during th e 



campaign and with previous studies of the X-ray activity of Mrk501 (IKataoka et al.ll200ll ). 
At this point, it is necessary to determine whether an emission region characterized by these 
values of R and 5 is optically thin to internal two-photon pair creation 77 — > e + e~ for the 
highest TeV energies observed during the campaign. We now affirm pair transparency due 
to insufficient density of soft target photons. 

Since Mrk 501 is a cosmologically local object, pair conversion in the EBL is not expected 
to prevent its multi-TeV photons from reaching the Earth, although the impact of this process 
is not negligible, as mentioned in £j5j Therefore, dealing with a nearby source allows us to 
focus mostly on the intrinsic absorption processes, rather than on the cosmological, EBL- 
related, attenuation of the 7-ray emission. Moreover, because of the absence (or weakness) 
of accretion-disk- related circumnuclear photon fields in BL Lac objects like Mrk 501, we 
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only need to consider photon-photon pair production involving photon fields internal to 
the jet emission site. The analysis is therefore simpler than in the case of FSRQs, where 
the attenuation of high-energy 7-ray fluxes is dominated by interactions with photon fields 
external to the jet — such as those provided by the broad line regions or tori — for which 
the exact spatial distribution is still under debate. 

Pair-creation optical depths can be estimated as follow s. Using the ^-function ap - 



proximation for the photon-photon annihilation cross-section ( jZdziarski fc Lightmanlll985l ). 
a^(e' Q , e' ) ~ 0.2 <7 T e' Q 5[e' Q — {2m 2 e c A / e'^)], the corresponding optical depth for a 7-ray photon 
with observed energy e 7 = 5 e' /(l + z) interacting with a jet-originating soft photon with 
observed energy 

£ o = 7— ^ ^ ,., , e N9 - 50 — -77 eV , (7) 



1 + z e 7 (l + zf V 10 / VTeV 
may be found as 

r 77 ~ ds de n' (e ) cr 77 (e ,£ 7 ) ~ 0.2 a T Re' n' (e ) , (8) 

J Jm e c 2 /eo 

where n' Q (e' Q ) is the differential comoving number density of soft photons. Noting that 
£'qTi' (£' ) = L' Q /47r R 2 c, where L' is the intrinsic monochromatic luminosity at photon energy 
e' , we obtain 



10 Rm 2 e c 5 5 5 ' VTeV/ V10" 11 erg/cm 2 /sy V10 17 cmy V 10 



where F = L /47rd 2 L is the observed monochromatic flux energy density as measured at 
the observed photon energy s . Thus, for 5TeV 7-rays and the model parameters discussed 
(implying the observed Eq = 15 eV flux of Mrk501 roughly F ~ 3.2 x 10~ n ergs^ 1 cm -2 ), 
one has r 77 (5TeV) ~ 0.005. Therefore, the values of R and 5 from our SSC model fit do 
not need to be adjusted to take into account the influence of spectral modifications due 
to pair attenuation. Note that such opacity effects, studied extensively in the context of 
7-ray bursts, generally yield a broken power law for the spectral form, with the pos ition 



and magnitude of the break fixed by the pair-production kinematics (e.g., iBaringi 120061 . and 
references therein). The broad-band continuum of Mrk501, and in particular its relatively 
flat spectrum VHE 7-ray segment, is inconsistent with such expected break. This deduction 
is in agreement with the above derived transparency of the emitting region for high energy 
7-ray photons. 

Next we evaluated the 'monoenergetic' comoving energy density of ultrarelativistic elec- 
trons for a given electron Lorentz factor, 

1 U' e ( 1 )= 1 2 m e c 2 n' e ( 1 ), (10) 
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and this is shown in FigurelTTI (solid black line). The total electron energy density is then 
U' e = J U'J^y) ~ 5 x 10 _4 ergcm~ 3 . As shown, most of the energy is stored in the lowest 
energy particles (7 m j n — 600). For comparison, the comoving energy density of the magnetic 
field and that of the synchrotron photons are plotted in the figure as well (horizontal solid 
red line and dotted blue line, respectively). These two quantities are approximately equal, 
namely U' B = B 2 /8tt ~ 0.9 x 10 _5 ergcm -3 and 



Ait R i 

U' sy n = -J^ I fu>,syndv' * 0.9 x lO^ergcm" 3 . (11) 



In Figured]] we also plot the comoving energy density of synchrotron photons which are 
inverse-Compton upscattered in the Thomson regime for a given electron Lorentz factor 7, 



AttR 



v(7) 



U'syn/Al) = J j'v',syn dv ' ( 12 ) 

(dashed blue line), where v'ki^i) = rn e c 2 /4ry h. Because of the well-known suppression of the 
inverse-Compton scattering rate in the Klein-Nishina regime, the scattering in the Thomson 
regime dominates the inverse-Compton energy losses^. Hence, one may conclude that even 
though the total energy density of the synchrotron photons in the jet rest frame is comparable 
to the comoving energy density of the magnetic field {U' syn ~ U' B ), the dominant radiative 
cooling for all the electrons is due to synchrotron emission, since U' syn , T < U B for any 7. 

The timescale for synchrotron cooling may be evaluated as 



4n-7^%-4f^)"day. (13) 



syn AcjtiU' \WJ 



Hence, t' rad ~ t' syn equals the dynamical timescale of the emitting region, t' dyn ~ R/c, for 
electron Lorentz factor 7 ~ 8 x 10 5 , i.e., close to the second electron break energy 7^2- 
Also the difference between the spectral indices below and above the break energy 7^2 
determined by our modeling, namely AS3/2 = S3 — S2 = 0.95, is very close to the 'classical' 
synchrotron cooling break As = 1 expected for a uniform emission region, as discussed in 
§6.11 This agreement, which justifies at some level the assumed homogeneity of the emission 
zone, was in fact the additional constraint imposed on the model to break the degeneracy 
between the main free parameters. Note that in such a case the first break in the electron 
energy distribution around electron energy ^br,i = 4 x 10 4 is related to the nature of the 
underlying particle acceleration process. We come back to this issue in §7.2[ 



The inverse-Compton cross-section goes a s a,; ff ~ <jt for v' < y 'v^'y)-, and roughly as <x 



cr T {v'/v'kn)' 1 Hf'/v'kn] for v ' > v'kn(i) (e.g- ICoppj fc Blandfordlll99(t 
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Another interesting result from our model fit comes from the evaluation of the mean 
energy of the electrons responsible for the observed non-thermal emission of Mrk501. In 
particular, the mean electron Lorentz factor is 

J <(7) «7 

This value, which is determined predominantly by the minimum electron energy jmin = 600 
and by the position of the first break in the electron energy distribution, is comparable to the 
proton-to-electron mass ratio m p /m e . In other words, the mean energy of ultrarelativistic 
electrons within the blazar emission zone of Mrk501 is comparable to the energy of non- 
relativistic/mildly-relativistic (cold) protons. This topic will be discussed further in §7.21 as 
well. 

The analysis presented allows us also to access the global energetics of the Mrk501 jet. 
In particular, with the given energy densities U' e and U' B , we evaluate the total kinetic powers 
of the jet stored in ultrarelativistic electrons and magnetic field as 

L e = ttR 2 cT 2 U' e ~ 10 44 erg s^ 1 , (15) 

and 

L B = ttR 2 cT 2 U' b ~ 2 x 10 42 erg s" 1 , (16) 

respectively. In the above expressions, we have assumed that the emission region analyzed 
occupies the whole cross-sectional area of the outflow, and that the jet propagates at suf- 
ficiently small viewing angle that the bulk Lorentz factor of the jet equals the jet Doppler 
factor emerging from our modeling, T = 5. These assumptions are justified in the frame- 
work of the one-zone homogeneous SSC scenario. Moreover, assuming one electron-pro ton 



pair per electron-positron pair within the emission region (see lCelotti fc Ghisellinill2008l ). or 
equivalently N p ~ N' e /3 where the total comoving number density of the jet leptons is 

we obtain the comoving energy density of the jet protons U' p = (j p ) m p c 2 N^/3, and hence 
the proton kinetic flux L p = ttR 2 cT 2 U p ~ 0.3 ( r y p ) 10 44 ergs _1 . This is comparable to the 
kinetic power carried out by the leptons for mean proton Lorentz factor (j p ) ~ 4 (see eq. [T5|) . 
It means that, within the dominant emission zone of Mrk501 (at least during non-flaring 
activity), ultrarelativistic electrons and mildly- relativistic protons, if comparable in number, 
are in approximate energy equipartition, and their energy dominates that of the jet magnetic 
field by two orders of magnitude. It is important to compare this result with the case of 
powerful blazars of the FSRQ type, for which the relatively low mean energy of the radiating 
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electrons, (7) <C 10 3 , assures dynamical do mination of cold pr otons even for a smaller proton 



content N'/N' e < 0.1 (see the discussion in lSikora et al.ll2009l and references therein). 



L p + Lb 

LEdd 



Assuming (7 p ) ~ 1 for simplicity, we find that the implied total jet power Lj = L e + 
1.4 x 10 44 ergs _1 constitutes only a small fraction of the Eddington luminosity 
47r GM B h Tripc /ar ^ (1.1-4.4) x 10 47 erg s" 1 for the Mrk 501 black hole mass M BH ~ 



(0.9 - 3.5) x 1O 9 M (IBarth et al.ll2002f ). In particular, our model implies Lj/L Edd ~ 10 



-3 



in Mrk 501. In this context , detailed investiga tion of the emission-line radiative output 
from the Mrk 501 nucleus by iBarth et al.l ( 120021 ) allowed for an estimate of the bolometric, 
accretion-related luminosity as L dis k — 2.4 x 10 43 ergs _1 , or Ldisk/LEdd ~ 10~ 4 . Such a 
relatively low luminosity is not surprising for BL Lac s, which are believed to accrete at 
low, sub-Eddington rates (e.g.. iGhisellini et al.l l2009bl ). For low-accretion-rate AGN (i.e., 
those for which L disk J ' L F/ m < 10~ 2 ) the expected radiative efficiency of the accretion disk is 
rjdisk = Ldi S k/L acc < 0.1 (INarayan fc Yilll994j ; iBlandford fc Begelmanl|l999l ). Therefore, the 
jet power estimated here for Mrk 501 is comfortably smaller than the available power of the 
accreting matter L acc . 



Finally, we note that the total emitted radiative power is 
L em ~ T 2 (L' n + L' ssc ) = AttR 2 cT 2 (U' + U' ssc ) 



10 43 erg s _1 , 



where U' syn is given in Equation [TT] and the comoving energy density of 7-ray photons, 
U' ssc , was evaluated in an analogous way as ~ 1.7 x 10~ 6 ergcm -3 . This implies that the 
jet/blazar radiative efficiency was at the level of a few percent (L em /Lj ~ 0.07) during the 
period covered by the multifrequency campaign. Such a relatively low radiative efficiency 
is a common characteristic of blazar jets in general, typic ally claimed to be at the level of 
1% — 10% fsee lCelotti k Ghisellinill2008Hsikora et al.ll2009h . On the other hand, the isotropic 



synchrotron and SSC luminosities of Mrk 501 corresponding to the observed average flux 
levels are L syn = 5 A L' syn ~ 10 45 ergs _1 and L ssc = 5 A L' SSC ~ 2 x 10 44 ergs _1 , respectively. 



7.2. Electron Energy Distribution 

The results of the SSC modeling presented in the previous sections indicate that the 
energy spectrum of freshly accelerated electrons within the blazar emission zone of Mrk 501 
is of the form oc E~ 2 - 2 between electron energy E e ^ min = ■y rn i n m e c 2 ~ 0.3 GeV and E £jbr = 
Ibr, i m e c 2 ~ 20 GeV, steepening to oc E~ 2 - 7 above 20 GeV, such that the mean electron en- 
ergy is (E e ) = (7) m e c 2 ~ 1 GeV. At this point, a natural question arises: is this electron 
distribution consistent with the particle spectrum expected for a diffusive shock accelera- 
tion process? Note in this context that the formation of a strong shock in the innermost 
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parts of Mrk501 might be expected around the location of the large bend (change in the 
position angle by 90°) observed in the outflow with in the first few parsecs (projected) from 
the core ( [Edwards fc Pinerl |2002| ; iPiner et al.l 120091 ) . This distance scale could possibly be 
reconciled with the expected distance of the blazar emission zone from the center for the 
model parameters discussed, r ~ R/0j ~ 0.5 pc, for jet opening angle ^ ~ 1/r C 1. 

In order to address this question, let us first discuss the minimum electron energy im- 
plied by the modeling, E e ^ min ~ 0.3 GeV. In principle, electrons with lower energies may 
be present within the emission zone, although their energy distribution has to be very flat 
(possibly even inverted), in order not to overproduce the synchrotron radio photons and to 
account for the measured Fermi-LAT 7-ray continuum. Therefore, the constrained minimum 
electron energy marks the injection threshold for the main acceleration mechanism, meaning 
that only electrons with energies larger than E e ^ min are picked up by this process to form 
the higher-energy (broken) power-law tail. Interestingly, the energy dissipation mechanisms 
operating at the shock fronts do introduce a particular characteristic (injection) energy scale, 
below which the particles are not freely able to cross the shock front. This energy scale de- 
pends on the shock microphysics, in particular on the thickness of the shock front. The shock 
thickness, in turn, is determined by the operating inertial length, or the diffusive mean free 
path of the radiating particles, or both. Such a scale depends critically on the constituents 
of the shocked plasma. For pure pair plasmas, only the electron thermal scale enters, and 
this sets E e ^ min ~ Tm e c 2 . In contrast, if there are approximately equal numbers of electrons 
and protons, the shock thickness can be relatively large. Diffusive shock acceleration can 
then operate on electrons only above a relatively high energy, establishing E e ^ min ~ eTm p c 2 . 
Here, e represents some efficiency of the equilibration in the shock layer between shocked 
thermal protons and their electron counterparts, perhaps result ing from electrostatic poten - 
tials induced by charge separation of species of different masses (IBaring fc Summerlinll2007l ). 
Our multifrequency model fits suggest that e ~ 0.025, providing an important blazar shock 
diagnostic. 

At e ven lower electron e nergies, other energization processes must play a dominant 
role (e.g., Hoshino et al.l 1992), resultin g in formation of very flat electron spectra (see the 
related discussion in lSikora et al.l 120091 ). Which of these energy dissipation mechanisms are 
the most relevant, as well as how flat the particle spectra could be, are subjects of ongoing 
debates. Different models and numerical simulations presented in the literature indicate 

from 



a wide possible range for the lowest-energy particle spectral indices (below E e 



1.0 — 1.5 down to s 



parameters (lAmato fc Arond 120061 ; ISironi fc Spitkovskyl 120091 ) . 



2, depending on the pa rticular shock conditions and 



All in all, we argue that the relatively high minimum energy of the radiating electrons 
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implied by the SSC modeling of the Mrk501 broadband spectrum and the overall character 
of the electron energy evolution in this source are consistent with a proton-mediated shock 
scenario. In addition, the fact that the mean energy of the ultrarelativistic electrons is of the 
order of the proton rest energy, (E e ) ~ m p c 2 , can be reconciled with such a model. Moreover, 
the constrained power-law slope of the electrons with energies E^ min < E e < E e ^ ri namely 
si = 2.2, seems to suggest a dominant role for diffusive shock acceleration above the injection 
energy E e>mini as this value of the spectral index is often claimed in the literature for par- 
ticles undergoing first-order Fermi acceleratio n at relativistic shock ( jBednarz fc Ostrowski 
19981 ; iKirk et al.l I200CH : lAchterberg et all l200lh . The caveat here is that this result for the 
'universal' particle spectral index holds only for particular conditions (namely, for ultrarel- 
ativistic shoc k velocities with highly turb ulent conditions downstream of the shock: see the 
discussion in lOstrowski fc Bednarzl |2002| ). whereas in general a variety of particle spectra 
may result from the relativistic first-order Fermi mechanism, depending on the local mag- 
netic field and turbulence parameters at the shock front, and the speed of the upstream flow 
( Kirk &: Heayenslll989uNiemiec fc Ostrowskill2004l ; iLemoine et al.ll2006t ISironi fc Spitkovsky 
20091 ; lBaringil2010l ). Nevertheless, the evidence for ultrarelativistic electrons with spectral 
index s\ = 2.2 in the jet of Mrk501 may be considered as an indication that the plasma 
conditions within the blazar emission zone allow for efficient diffusive shock acceleration (at 
least in this source), as described by the simplest asymptotic test-particle models, though 
only in a relatively narrow electron energy range. 

If the relativistic shock acceleration plays a dominant role in the blazar Mrk 501 (as 
argued above), the observations and the model results impose important constraints on this 
mechanism, many aspects of which are still not well understood. Firstly, this process must 
be very efficient in the sense that all the electrons pre- accelerated/preheated to the energy 
E e ^min ~ 0.3 GeV are picked up by the acceleration process so that a single electron compo- 
nent is formed above the injection threshold and there is no Maxwellian-like population of 
particles around E e ^ min outnumbering the higher energy ones from the power-law taiQ The 
second constraint is due to the presence of the spectral break As 2 /i = $2 — si — 0.5 around 
electron energies E e , b r ~ 20 GeV. As discussed in the previous section, this break cannot be 
simply a result of cooling or internal pair-attenuation effects, and hence must be accounted 
for by the acceleration mechanism. The discussion regarding the origin of this break - 
which may reflect variations in the global field orientation or turbulence levels sampled by 



24 Note that due to low accretion rates and thus low luminosities of the accretion disks in BL Lacs, the 
number of non-relativistic/mildly- relativistic electrons (Lorentz factors 7 ~ 1) cannot be constrained by 
analyzing 'b ulk-Compton' spectral fe atures in the observed SED of Mrk 501, in contrast to the situation in 
FSRQs Csee lSikora fc Madeiskillioooh . 
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particles of different energy — is beyond the scope of this paper. However, the presence 
of high-energy breaks in the electron energy distribution (intrinsic to the particle spectrum 
rather than forming due to cooling or absorption effects) seems to be a common property 
of relativistic j et sources observed by Fe rmi-LAT, such as the FSRQ objects 3C 454.3 and 
AO 0235+164 (jAbdo et al.ll2009bl . boiOah . 



7.3. Variability 

In §3]-§U we reported on the 7-ray flux and spectral variability of Mrk501 as measured 
by the Fermi-LAT instrument during the first 16 months of operation. In this section, we 
discuss whether the observed spectral evolution can be accounted for by our simple one-zone 
SSC model. 

Figured^] (top panel) presents in more detail the GeV-TeV 7-ray spectrum of Mrk501, 
together with the decomposition of the SSC model continuum. Here the contributions of 
different segments of the electron energy distribution are indicated by different colors. As 
shown, the low-energy electrons, 7 m j n < 7 < 7^ 1, which emit synchrotron photons up to 
~ 10 15 Hz, dominate the production of 7-rays up to a few GeV (red line). The contribution 
of higher-energy electrons with Lorentz factors 7^ 1 < 7 < 7^ 2 is pronounced within the 
observed synchrotron range 10 15 — 10 18 Hz, and at 7-ray energies from a few GeV up to ~ TeV 
(green line). Finally, the highest energy tail of the electron energy distribution, 7 > 7^2, 
responsible for the observed hard- X-ray synchrotron continuum (> 2keV) in the fast cooling 
regime, generates the bulk of 7-rays with observed energies > TeV (blue line). Interestingly, 
even though any sharp breaks in the underlying electron energy distribution are 'smeared 
out' into a smoothly curved spectral continuum due to the nature of the synchrotron self- 
Compton emission, the average data set does support the presence of distinct low-energy 
and high-energy segments in the electron spectrum. 

It therefore seems reasonable to argue that the spectral variability of Mrk501 observed 
by Fermi-LAT may be explained by postulating that the low-energy segment of the electron 
energy distribution (7 < 7^,1) is characterized by only small flux variations, while the 
high-energy electron tail (7 > 7& r ,i) varies more substantially. In such a scenario, some 
correlation might be expected between the fluxes in the UV-to-soft-X-ray photon energies 
from the synchrotron bump and the GeV-TeV fluxes from the inverse- Compton bump. This 
expectation is not inconsistent with the fact that we do not see any obvious relation between 
the ASM/BAT fluxes and the LAT (> 2 GeV) fluxes (see Figure|2]), because the electrons 
producing X-ray synchrotron photons above 2 — 3 keV contribute to the SSC emission mostly 
at the highest photon energies in the TeV range (see Figure [T2"j) . This issue will be studied 
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in detail (on timescales of 1 month down to 5 days) in a forthcoming publication using the 
data from this multifrequency campaign. 

The X-ray/TeV connection has been established in the past for many BL Lacs (and for 
Mrk501 in particular). However, the exact character of the correlation between the X-ray 
and TeV fluxes is known to vary from object to object, a n d from epoch to epoch in a give n 



source, as w i dely d i scussed in, e . g.. iKrawczynski et al.l (120041 ); iBlazejowski et al.l ( 120051 ): 



Gliozzi et al.l ( 120061 ); iFossati et al.l ( 120081 ). Note that the data analyzed in those papers were 
obtained mostly during periods of high activity. Consequently, the conclusions presented 
were somewhat biased towards flaring activity, and hence they might not apply to the typical 
(average) behaviour of the source, which is the main focus of this paper. Moreover, the data 
set from our campaign includes UV fluxes, soft-X-ray fluxes (down to 0.3 keV; Swift), and 7- 
ray fluxes spanning a very wide photon energy range (0.1 GeV— 10 TeV; Fermi-LAT combined 
with MAGIC and VERITAS). This unique data set allows us to evaluate the multifrequency 
variability and correlations for Mrk501 over an unprecedented range of photon energies. 

Considering only the data set presented in this paper, we note that by steepening the 
high-energy electron continuum above the intrinsic break energy 7^1 (and only slightly 
adjusting the other model parameters), one can effectively remove photons above 10 GeV in 
the SSC component, leaving a relatively steep spectrum below 10 GeV, similar to the one 
observed by Fermi-LAT during the time interval MJD 54862-54892 (see §4]). Such a change 
should be accompanied by a decrease in the UV-to-soft-X-ray synchrotron fluxes by a factor 
of a few, but the data available during that time interval are not sufficient to detect this 
effectP^I. This statement is further justified by the bottom panel in Figure[T2l where the 
contributions of the different segments of electrons Comptonizing the different segments of 
the synchrotron bump to the average 7-ray emission of Mrk501 are shown. Note that the 
lowest-energy electron population (7 < 7& r> i) inverse- Compton upscattering only synchrotron 
photons emitted by the same population [y < z/ br> 1 ~ 10 15 Hz; solid red line) may account 
for the bulk of the observed steep-spectrum 7-ray emission. 

Another important conclusion from this figure is that Comptonization of the highest- 
energy synchrotron photons [y > f& r ,2 ~ 6 x 10 17 Hz) by electrons with arbitrary energies 
produces only a negligible contribution to the average 7-ray flux of Mrk501 due to Klein- 
Nishina suppression. Thus, the model presented here explains in a natural way the fact 
that the X-ray and TeV fluxes of TeV-emitting BL Lacs are rarely correlated according to 



25 The 4. 5-month multifrequency campaign started 13 days after the end of the 30-day time interval MJD 
54862-54892. Therefore, for this epoch, the only additional multifrequency data are from RXTE-ASM and 
Swift-BAT, which have only moderate ability to detect Mrk501 on short timescales. 
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the simple scaling i*r e v oc F^ eV which would be expected from the class of SSC models in 
which the highest-energy electrons u pscatter (in the Thomson regime) their own synchrotron 
photons to the TeV band (see, e.g., Gliozzi et al. 2006 ; Fossati et al.l 2008 ). In addition, it 
opens a possibility for accommodating short-timescale variability (t var < 4d) at the high- 
est synchrotron and inverse-Compton frequencies (hard X-rays and TeV photon energies, 
respectively). The reason for this is that, in the model considered here, these high-energy 
tails of the two spectral components are produced by the highest-energy electrons which 
are deep in the strong cooling regime (i.e., for which t' rad <C R/c), and thus the corre- 
sponding flux changes may occur on timescales sh orter than R/c5 (see in this context, e.g., 



Chiaberge fc Ghisellinilll999l : Kataoka et al.ll2000h . 



A more in-depth analysis of the multifrequency data set (including correlation studies 
of the variability in different frequency ranges) will be presented in a forthcoming paper. 
The epoch of enhanced 7-ray activity of Mrk501 (MJD 54952-54982; see §4|) may be more 
difficult to explain in the framework of the one-zone SSC model, because a relatively flat 
Fermi-LAT spectrum above 10 GeV, together with an increased TeV flux as measured by the 
VERITAS and Whipple 10 m telescopes around this time, may not be easy to reproduce with 
a set of model parameters similar to that discussed in previous sections. This is mostly due to 
Klein-Nishina effects, which tend to steepen the high-energy tail of the SSC component, thus 
precluding the formation of a flat power law extending beyond the observed TeV energies. 
Hence, detailed modeling and data analysis will be needed to determine whether the enhanced 
VHE 7-ray activity period can be accommodated within a one-zone SSC model, or whether 
it will require a multi-zone approach. 



8. Conclusions 



We have presented a study of the 7-ray activity of Mrk501 as measured by the LAT 
instrument on board the Fermi satellite during its first 16 months of operation, from 2008 
August 5 (MJD 54683) to 2009 November 27 (MJD 55162). Because of the large leap in 
capabilities of LAT with respect to its predecessor, EGRET, this is the most extensive study 
to date of the 7-ray activity of this object at GeV- TeV photon energies. The Fermi-LAT 
spectrum (fitted with a single power-law function) was evaluated for 30-day time intervals. 
The average photon flux above 0.3 GeV was found to be (2.15±0.11) x 10~ 8 phcm~ 2 s _1 , and 
the average photon index 1.78 ± 0.03. We observed only relatively mild (factor less than 2) 
7-ray flux variations, but we detected remarkable spectral variability. In particular, during 
the four consecutive 30-day intervals of the "enhanced 7-ray flux" (MJD 54862-54982), the 
photon index changed from 2.51 ± 0.20 (for the first interval) down to 1.63 ± 0.09 (for the 
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Table 2. Parameters of the blazar emission zone in Mrk501. 



Parameter 


Main SSC fit considered 




Exemplary SSC fit 


Magnetic field 


B = 


0.015 G 


B = 


0.03 G 


Emission region size 


R = 


1.3 x 10 17 cm 


R = 


0.2 x 10 17 cm 


Jet Doppler and bulk Lorcntz factors 


r = s = 12 


r = , 


5 = 22 


Equipartition parameter 


Ve = 


U'JU' B = 56 




U'JUL = 130 

e/ B 


Minimum electron energy 


"tmin 


= 600 


' ' 


= 300 


Intrinsic electron break energy 


Ibr, 1 


= 4 x 10 4 


76r 1 


= 3 x 10 4 


Cooling electron break energy 


Ibr, 2 


= 9 x 10 5 


Ibr, 2 


= 5 x 10 5 


Maximum electron energy 




= 1.5 x 10 7 


^maa 


: =0.3x 10 7 


Low-energy electron index 


51 = 


2.2 


Si = 


2.2 


High-energy electron index 


S'2 = 


2.7 


S'2 = 


2.7 


Electron index above the cooling break 


S3 = 


3.65 


S3 = 


3.5 


Mean electron energy 


<7> = 


i 2400 


< 7 >, 


i 1200 


Main variability timcscale 


tvar 


~ 4 day 


tvar 


~ 0.35 day 


Comoving electron energy density 




0.5 X 10~ 3 ergcm~ 3 




: 4.6 X 10~ 3 crgcm~ 3 


Comoving magnetic field energy density 


U' B , 


; 0.9 x 10~ 5 erg cm -3 


U' B , 


i 3.6 X 10~ 5 crgcm~ 3 


Comoving energy density of synchrotron photons 


U' 

^ syn 


~ 0.9 X 10 _s erg cm" 3 


U' 

^ syn 


~ 3.1 X 10~ 5 crgcm~ 3 


Comoving electron number density 


N' e - 


; 0.3 cm- 3 


N' e - 


; 4.6 cm- 3 


Luminosity of the host galaxy 


Lstai 


. ~ 3 X 10 44 ergs" 1 


Lsta', 


- ~ 3 X lO^crgs- 1 


Jet power carried by electrons 


L e ~ 


1.1 X lO^ergs- 1 


L e ~ 


: 0.85 X 10 44 ergs" 1 


Jet power carried by magnetic field 


L B - 


t 2 X 10 42 ergs" 1 


L B - 


t 0.65 X 10 42 ergs" 1 


Jet power carried by protons™ 


Lp — 


3 X 10 43 ergs" 1 


Lp 2: 


: 4.2 X 10 43 ergs" 1 


Total jet kinetic power 


Lj- 


1.4 X lO^crgs- 1 


Lj- 


1.3 X 10 44 ergs" 1 


Total emitted power 


Lp/m 


~ 9.7 X 10 42 ergs" 1 


L em 


~ 2.7 X 10 42 ergs" 1 


Isotropic synchrotron luminosity 


Lgyn 


~ lO^'crgs- 1 


Lsyn 


~ 10 45 ergs" 1 


Isotropic SSC luminosity 


L S sc 


~ 2 X 10 44 ergs" 1 


L S sc 


~ 2 X lO^crgs- 1 



Note. — a Assuming one electron-proton pair per electron-positron pair, and mean proton Lorentz factor {■jp) ~ 1. 
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fourth one). During the whole period of 16 months, the hardest spectral index within the 
LAT energy range was 1.52 ± 0.14, and the softest one was 2.51 ± 0.20. Interestingly, this 
outstanding (and quite unexpected) variation in the slope of the GeV continuum did not 
correlate with the observed flux variations at energies above 0.3 GeV. 

We compared the 7-ray activity measured by LAT in two different energy ranges (0.2 — 
2 GeV and > 2 GeV) with the X-ray activity recorded by the all-sky instruments RXTE- 
ASM (2 — lOkeV) and Swift-BAT (15 — 50keV). We found no significant difference in the 
amplitude of the variability between X-rays and 7-rays, and no clear relation between the 
X-ray and 7-ray flux changes. We note, however, that the limited sensitivity of the ASM 
and (particularly) the BAT instruments to detect Mrk501 in a 30-day time interval, together 
with the relatively stable X-ray emission of Mrk501 during the observations, precludes any 
detailed X-ray/7-ray variability or correlation analysis. 

In this paper we also presented the first results from a 4.5-month multifrequency cam- 
paign on Mrk501, which lasted from 2009 March 15 (MJD 54905) to 2009 August 1 (MJD 
55044). During this period, the source was systematically observed with different instru- 
ments covering an extremely broad segment of the electromagnetic spectrum, from radio 
frequencies up to TeV photon energies. In this manuscript, we have focused on the average 
SED emerging from the campaign. Further studies on the multifrequency variability and 
correlations will be covered in a forthcoming publication. 

We have modeled the average broadband spectrum of Mrk501 (from radio to TeV) 
in the framework of the standard one-zone synchrotron self-Compton model, obtaining a 
satisfactory fit to the experimental data. We found that the dominant emission region in 
this source can be characterized by a size of R ~ 10 3 r 9 , where r g ~ 1.5 x 10 14 cm is the 
gravitational radius of the black hole (Mbh — 1O 9 M ) hosted by Mrk501. The intrinsic (i.e., 
not affected by cooling or absorption effects) energy distribution of the radiating electrons 
required to fit the data was found to be of a broken power-law form in the energy range 
0.3 GeV— 10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 
E B: br ~ 20 GeV, respectively. In addition, the model parameters imply that all the electrons 
cool predominantly via synchrotron emission, forming a cooling break at 0.5 TeV. We argue 
that the particular form of the electron energy distribution emerging from our modeling is 
consistent with the scenario in which the bulk of the energy dissipation within the dominant 
emission zone of Mrk501 is related to relativistic, proton-mediated shock waves. The low- 
energy segment of the electron energy distribution [E e < E e ^ r .) formed thereby, which 
dominates the production of 7-rays observed below a few GeV, seems to be characterized 
by low and relatively slow variability. On the other hand, the high-energy electron tail 
(E e > E ejbr ), responsible for the bulk of the 7-rays detected above a few GeV, may be 
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characterized by more significant variability. 

Finally, we found that ultrarelativistic electrons and mildly-relativistic protons within 
the blazar zone of Mrk 501, if comparable in number, are in approximate energy equipartition, 
with their energy dominating the energy in the jet magnetic field by about two orders of 
magnitude. The model fit implies also that the total jet power, Lj ~ 10 44 ergs _1 , constitutes 
only a small fraction of the Eddington luminosity, Lj/LEdd ~ 10~ 3 , but is an order of 
magnitude larger than the bolometric, accretion-related luminosity of the central engine, 
Lj/Ldi S k ~ 10. Finally, we estimated the radiative efficiency of the Mrk 501 jet to be at the 
level of a few percent, L em /Lj < 0.1, where L em is the total emitted power of the blazar. 
The results from this study could perhaps be extended to all HSP BL Lacs. 
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Fig. 7. — Spectral energy distribution for Mrk501 from Fermi-LAT for several time intervals 
of interest. The panel (a) shows the SED for the time period before the multifrequency 
campaign (MJD 54683-54901), the panel (b) for the time interval corresponding to the 
multifrequency campaign (MJD 54905-55044) and the panel (c) for the period after the 
campaign (MJD 55044-55162). In all panels, the black line depicts the result of the unbinned 
likelihood power-law fit, the red contours denote the 68% uncertainty of the power-law fit 
and blue arrows denote upper limits at 95% confidence level, which were computed for the 
differential energy ranges with a signal of TS < 4 or less than two photons predicted by 
the analysis model for Mrk501. The legend reports the results from the unbinned likelihood 
power-law fit in the energy range 0.3 — 400 GeV. 
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Fig. 8. — Spectral energy distribution for Mrk501 averaged over all observations taken during 
the multifrequency campaign performed between 2009 March 15 (MJD 54905) and 2009 
August 1 (MJD 55044). The legend reports the correspondence between the instruments 
and the measured fluxes. Further details about the instruments are given in §5.11 The 
optical and X-ray data have been corrected for Galactic extinction, but the host galaxy 
(which is clearly visible at the IR/optical frequencies) has not been subtracted. The TeV 
data from MAGIC and VERITAS have been corrected for the absorpti on in the extragalactic 
background light using the model reported in iFranceschini et al.l (120081 ) . The VERITAS data 
from the time interval MJD 54952.9-54955.9 were removed from the data set used to compute 
the average spectrum, and are depicted separately in the SED plot (in green diamonds). See 
text for further details. 
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Fig. 9. — Top panel: Enlargement of the 7-ray energy range from Figured Bottom 
panel: Same SED as in the top panel, but with the Fermi-LAT data from the multifrequency 
campaign split in two data sets: MJD 54952-54982 (open blue squares) and the rest (filled 
blue circles). 
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Fig. 10. — The SSC model fits to the broadband emission spectrum of Mrk501, averaged 
over all the observations made during the multifrequency campaign performed between 2009 
March 15 (MJD 54905) and 2009 August 1 (MJD 55044). The red bow-tie in the figure 
corresponds to the 68% containment of the power-law fit to the average Fermi-LAT spectrum 
(photon index 1.74 ± 0.05). The dotted black curve denotes the fit to th e starlight emissio n 
of the host galaxy assuming a template of a luminous elliptical as given in lSilva et al.l (119981 ). 
The details of the model are given in £j6l The black curves correspond to the main set of the 
model parameters considered (variability timescale t var ~ 4 days), while the red dot-dashed 
curves to the alternative set of the model parameters with the emission region size decreased 
by an order of magnitude (t var ~ 0.35 day). See text for further discussion. 
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Fig. 11. — The jet comoving energy density of ultrarelativistic electrons per logarithmic 
energy bin, 7^(7), as a function of the electron Lorentz factor 7 (solid black curve). For 
comparison, the comoving energy density of the magnetic field U' B (solid red line) and syn- 
chrotron photons U' syn (dotted blue line) are shown. The dashed blue curve denotes the 
comoving energy density of synchrotron photons which are inverse-Compton upscattered in 
the Thomson regime, U' , T , for a given electron Lorentz factor 7 (see equation [T2j) . 
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Fig. 12. — The decomposition of the SSC continuum for Mrk501. The data points are the 
same as in the bottom panel of Figure IHJ The SSC fit to the average spectrum is denoted 
by the solid black curve. Top: Contributions of the different segments of the electron 
spectrum Comptonizing the whole synchrotron continuum (red curve: ■y m i n < 7 < jbr,i] 
green curve: 7^1 < 7 < 7^,2; blue curve: 7^2 < 7). Bottom: Contributions of the 
different segments of the electron spectrum (as in the top panel) Comptonizing different 
segments of the synchrotron continuum (solid curves: v < Vbr,i — 10 15 Hz; dashed curves: 
Vbr, 1 < v < Vbr,2 — 6 x 10 17 Hz; curves for v > Vbr,2 do not appear in the plot because the 
corresponding flux levels are all less than 10 -13 erg cm -1 s^ 1 ). 



